NBSIR  80-2157 

?30713 

A Computer  Program  for  Analysis 
of  Pressurized  Stairwells  and 
Pressurized  Elevator  Shafts 


John  H.  Klote 


Center  for  Fire  Research 
National  Engineering  Laboratory 
U.S.  Department  of  Commerce 
National  Bureau  of  Standards 
Washington,  DC  20234 


January  1981 


Final  Report 


f-QG 

100 
I.U56 
80-2157 
1981 


U.S.  DEPARTMENT  OF  COMMERCE 
NATIONAL  BUREAU  OF  STANDARDS 


NBSIR  80-2157 

» * > 

A COMPUTER  PROGRAM  FOR  ANALYSIS 
OF  PRESSURIZED  STAIRWELLS  AND 
PRESSURIZED  ELEVATOR  SHAFTS 


National  Bureau  of  Standards 
Library,  E-01  Admin.  Bldg, 
FEB  2 7 1981 

f\Oi  &Ct-  - 

Clooo 

S'o 

fQSt 


John  H.  Klote 


Center  for  Fire  Research 
National  Engineering  Laboratory 
U.S.  Department  of  Commerce 
National  Bureau  of  Standards 
Washington,  DC  20234 


January  1981 
Final  Report 


Sponsored  in  part  by: 

The  Department  of  Health  and  Human  Services 
Washington,  DC  20201 


U.S.  DEPARTMENT  OF  COMMERCE,  Philip  M.  Klutznick,  Secretary 
Jordan  J.  Baruch,  Assistant  Secretary  for  Productivity.  Technology,  and  Innovation 
NATIONAL  BUREAU  OF  STANDARDS,  Ernest  Ambler.  Director 


PREFACE 


This  report  is  an  interim  product  of  a joint  effort  of  the  Department 
of  Health  and  Human  Services  and  the  National  Bureau  of  Standards  (NBS) , 
Center  for  Fire  Research.  The  program  is  a multi-year  activity  initiated  in 
1975.  It  consists  of  projects  in  the  areas  of:  decision  analysis,  fire  and 

smoke  detection,  smoke  movement  and  control,  automatic  extinguishment,  and 
behavior  of  institutional  populations  in  fire  situations. 

This  report  describes  a computer  program  which  analyzes  pressurized 
stairwells  and  pressurized  elevators.  The  program  was  initially  intended  as 
a research  tool  to  investigate  the  feasibility  of  specific  systems.  How- 
ever, this  program  may  be  of  interest  to  design  engineers  responsible  for 
pressurized  stairwells  or  pressurized  elevators. 
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A COMPUTER  PROGRAM  FOR  ANALYSIS  OF  PRESSURIZED 
STAIRWELLS  AND  PRESSURIZED  ELEVATOR  SHAFTS 

John  H.  Klote 
Abstract 

Pressurized  stairwells  and  pressurized  elevators 
can  be  used  as  a means  of  providing  a smoke  free  exit 
route  during  fire  situations.  This  paper  describes  a 
computer  program  which  analyzes  systems  intended  to 
pressurize  stairwells  or  elevator  shafts.  The  basic 
assumptions  and  limitations  of  the  program  are  also 
discussed.  The  appendices  contain  a program  listing 
and  examples. 

Key  words:  Air  movement;  computer  programs;  egress; 

elevator  shafts;  escape  means;  modeling;  pressurization; 
simulation;  smoke  control;  stairwells. 


1.  INTRODUCTION 

Pressurized  stairwells  are  being  used  as  a means  of  providing  smoke 
free  exit  routes  during  fire  situations.  A discussion  of  several  designs 
for  pressurized  stairwell  systems  is  provided  by  Benjamin  and  Klote  [1]  . 

In  addition,  the  concept  of  pressurized  elevator  shafts  has  gained  consider- 
able interest  as  a means  of  fire  escape  for  the  handicapped.  The  purpose  of 
the  program  described  in  this  paper  is  to  calculate  the  pressure  differentials 
produced  by  such  systems.  Because  a stairwell  or  an  elevator  shaft  is 
connected  to  a building,  the  program  calculates  air  flows  and  pressure 
differentials  throughout  the  building  in  order  to  obtain  the  differential 
pressures  across  the  stairwell  or  elevator. 

A number  of  computer  programs  have  been  developed  which  are  applicable 
to  smoke  control.  Some  of  these  programs  calculate  steady  state  air  flow 
and  pressures  throughout  a building  [2,3].  Other  programs  go  beyond  this  to 
calculate  smoke  concentrations  throughout  a building  that  would  be  produced 
in  the  event  of  a real  fire  [4-9] . In  general,  most  of  these  programs  are 
capable  of  analyzing  pressurized  stairwells  and  pressurized  elevators. 

However,  the  program  described  in  this  paper  has  been  specifically  written 
for  analysis  of  pressurized  stairwells  and  elevators.  The  data  input  has 


^"Numbers  in  brackets  refer  to  the  literature  references  listed  at  the  end  of 
this  paper. 
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been  designed  to  minimize  the  quantity  of  required  data  and  still  maintain  a 
high  level  of  generality  in  the  model.  The  output  consists  of  the  differen- 
tial pressures  across  all  of  the  building  shafts  in  addition  to  the  flows 
and  pressures  throughout  the  building.  In  addition,  the  assumptions  and 
limitations  of  the  program  discussed  in  the  next  section  were  chosen  specif- 
ically to  allow  modeling  of  pressurized  stairwells  and  pressurized  elevators. 

This  program  was  intended  primarily  as  a research  tool  to  investigate 
the  feasibility  of  specific  systems  and  to  determine  the  interaction  between 
these  systems  and  the  rest  of  the  building.  This  program  has  already  been 
used  to  analyze  pressurized  stairwells  without  vestibules  and  to  evaluate 
factors  which  affect  the  performance  of  these  systems  [10] . It  is  also 
possible  that  the  program  may  be  used  as  a design  tool  directly  or  to  generate 
quantities  of  design  data.  This  paper  is  not  intended  to  be  a design  guide 
for  either  pressurized  stairwells  or  pressurized  elevators.  The  state-of- 
the-art  of  these  systems  is  still  developing  and  designers  of  these  systems 
should  seek  the  most  current  data  available. 

2.  PROGRAM  CONCEPT 

In  this  computer  program  a building  is  represented  by  a network  of 
spaces  or  nodes  each  at  a specific  pressure  and  temperature.  The  stairwells 
and  other  shafts  are  modeled  by  a vertical  series  of  spaces,  one  for  each 
floor.  Air  flows  through  leakage  paths  from  regions  of  high  pressure  to 
regions  of  low  pressure.  These  leakage  paths  are  doors  and  windows  which 
may  be  opened  or  closed.  Leakage  can  also  occur  through  partitions,  floors, 
exterior  walls  and  roofs.  The  air  flow  through  a leakage  path  is  a function 
of  the  differential  pressure  across  the  leakage  path. 

In  this  computer  model  air  from  outside  the  building  can  be  introduced 
by  a pressurization  system  into  any  level  of  a shaft  or  even  into  other 
building  spaces.  This  allows  simulation  of  stairwell  pressurization,  eleva- 
tor shaft  pressurization,  stairwell  vestibule  pressurization  and  even  eleva- 
tor lobby  pressurization.  The  pressures  throughout  the  building  and  flow 
rates  through  all  the  flow  paths  are  obtained  by  solving  the  air  flow  net- 
work including  the  driving  forces  such  as  the  wind,  the  pressurization 
system  or  an  inside  outside  temperature  difference. 
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3.  ASSUMPTIONS  AND  LIMITATIONS 


1.  Each  space  is  considered  to  be  at  one  specific  pressure  and  one  specific 
temperature . 

2.  The  flows  and  leakage  paths  are  assumed  to  occur  at  mid-height  of  each 
level . 

3.  The  net  air  supplied  by  the  air  handling  system  or  by  the  pressurization 
system  is  assumed  to  be  constant  and  independent  of  building  pressure. 

4.  The  outside  air  temperature  is  assumed  to  be  constant. 

5.  The  barometer  pressure  at  ground  level  is  assumed  to  be  standard 
atmospheric  pressure  (101325  P ) . The  results  of  the  program  are  not 

3. 

very  sensitive  to  changes  in  atmospheric  pressure.  For  altitudes 
considerably  different  from  sea  level  the  more  accurate  value  can  be 
substituted  by  changing  an  assign  statement  in  the  subroutine  INPUT. 


4.  EQUATIONS 


A.  Flow  equation 


F = C (AP) x 


(3.1) 


where : 


F = mass  flow  rate 


C = flow  coefficient 


X = flow  exponent 


AP  = differential  pressure  across  flow  path 


The  following  equation  has  been  used  by  Sander  and  Tamura  [3]  for 

2 

corrections  for  flows  not  at  standard  conditions  . 


C - Co  <£-> 


X 


1 — 2 X 


o 


where : 


C = flow  coefficient  at  standard  conditions 


o 


p = density  of  entering  air 
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This  relationship  was  originally  derived  by  Tamura  and  Wilson  [11] . 
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PQ  = density  of  air  at  standard  conditions 

U = dynamic  viscosity  of  entering  air 

UQ  = dynamic  viscosity  of  air  at  standard  conditions 

This  author  has  observed  that  the  dynamic  viscosity  is  proportional  to 
absolute  temperature  to  the  3/4  power  in  the  range  of  -40°C  to  60°C  (-40°F 
to  140°F)  with  a maximum  error  of  1.2  percent.  Using  this  and  the  fact  that 
density  is  inversely  proportional  to  absolute  temperature  the  above  relation 
reduces  to  a function  of  temperature  only. 


where : 

T = absolute  temperature  of  entry  air 
T = standard  absolute  temperature 


B . Mass  Balance  Equations 


3 . 

For  building  compartment  l 


N 


N 


E Fd,j) + E 


j=i 


k=l 


F , . . . + F - ...  =0 

o(i,k)  f(i) 


(3.3) 


and  for  shafts 


i=N. 


N 


N 


S F(i,j)  +]C  Fo  (i,k)  + Ff(i] 
j=l 


k=l 


= 0 


where : 


(3.4) 


F (i / j ) 


mass  flow  rate  from  space  j to  space  i.  For  building 
compartments  this  flow  can  be  either  horizontal  or  vertical, 
however  for  shafts  this  flow  can  only  be  horizontal. 


o (i ,k) 


mass  flow  rate  from  direction  k outside  of  the  building  to 
space  i. 


3 

In  this  paper  the  term  building  compartment  refers  to  a space  in  a building 
other  than  in  a shaft. 
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Ff ( i)  = net  mass  flow  rate  of  air  due  to  the  air  handling  system  or 
due  to  a pressurization  system. 

Nc  = number  of  building  spaces  connected  to  space  i. 

N = number  of  connections  to  the  outside  from  space  i. 

o 

is  the  space  number  at  bottom  level  of  shaft  and  the  spaces  in  the 

shaft  are  numbered  consecutively  up  to  ^ which  is  the  space  number  at 

the  top  of  the  shaft. 

C.  Shaft  Pressures 

The  following  relation  is  used  to  calculate  the  gage  pressure,  P(i)'  at 
level  i of  a shaft  in  terms  of  P^_^  at  level  i - 1. 


P 


(i) 


P 


(i-D 


P 

z 


- P, 


(3.5) 


where : 

P = hydrostatic  pressure  difference 
P^  = pressure  loss  due  to  friction 

The  following  equation  is  used  to  calculate  the  hydrostatic  pressure. 


= 2Z 

RT 


hd)  " hd 


(3.6) 


where : 


h ^ = height  of  point  i 

h^_u  = height  of  point  i-1 

g = gravitational  acceleration 

R = gas  constant 

- _ T(i)  + T (i-1) 


_ P(i)  + P (i-1) 


+ P, 


2 b 

P^  is  a constant  used  to  convert  an  average  gage  pressure  to  the  average 
absolute  pressure  P. 


The  following  equation  is  used  to  calculate  the  pressure  loss  due  to  friction, 


5 


(3.7) 


where : 

= upward  flow  from  i-1  to  i in  shaft 

C = shaft  flow  coefficient 
s 

S = sign  of  F 
^ u 

D . Outside  Pressures 

Outside  pressures  can  either  be  entered  by  the  user  or  can  be  calculated  by 
the  following  method. 


P / ' \ ~ P-U  / • \ + C P , . » 
0(1)  h(i)  w v(i) 

where : 


(3.8) 


PQ^  = outside  gage  pressure  at  height  h(i)  above  absolute  pressure  at 
ground  level. 

Ph(i)  = hYdrostatic  Pressure  difference  between  h(i)  and  ground  level 

Pv^j  = velocity  pressure  due  to  the  wind  at  height  h(i) 

C = pressure  coefficient 

w 

Because  the  outside  temperature  is  constant 


where : 

Patm  = abs°lute  barometric  pressure  at  ground  level 
TQut  = outside  absolute  temperature 


When  the  outside  pressures  are  calculated  by  the  computer  the  wind 
velocities  are  assumed  to  be  described  by  the  power  law. 

V = V {£-)n 

o h 

o 

where : 

V = wind  velocity  at  height  h 

o o 

n = wind  exponent 
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This  relationship  has  been  extensively  used  to  describe  the  boundary- 
layer  velocity  profile  of  the  wind  near  the  surface  of  the  earth.  This 
equation  assumes  that  the  terrain  surrounding  the  building  is  homogeneous. 
That  is,  that  there  are  no  large  obstructions  near  the  building  which  could 
produce  local  wind  effects.  A value  of  0.16  for  the  wind  exponent  is  appro 
priate  for  flat  terrain.  The  wind  exponent  increases  with  rougher  terrain. 
For  very  rough  terrain  such  as  urban  areas  a value  of  0.40  would  be 
appropriate . 

The  equation  for  the  velocity  pressure  at  height  h(i)  is  obtained  by 
substituting  the  velocity  from  the  power  law  into  the  usual  relation  for 
velocity  pressure  (Pv  = j pV2) . 


where  p is  the  outside  air  density. 

The  pressure  coefficients  are  in  the  range  of  -0.8  to  0.8  where  positive 
values  are  for  windward  walls  and  negative  values  are  for  leeward  walls. 

The  pressure  coefficient  depends  upon  building  geometry  and  varies  locally 
over  the  wall  surface.  Numerical  values  for  and  n as  well  as  practical 
engineering  information  is  available  from  a number  of  sources  [12-15]  . 


5.  PROGRAM  DESCRIPTION 

This  program  is  written  in  the  computer  language  FORTRAN  V specifically 
for  use  on  the  UNIVAC  1108.  A program  listing  is  provided  in  appendix  D. 

The  following  is  a detailed  description  of  the  main  program  and  the  major 
subroutines . 


5.1  Main  Program 

This  program  was  written  to  handle  four  different  run  types  as  follows: 

1.  RUN:  In  this  run  type  subroutine  INPUT  is  used  to  read  data  and  then 
the  pressures  and  flows  are  calculated  and  printed.  After  this,  the 
subroutine  INPUT2  is  used  to  read  modifications  to  existing  data  or  to 
determine  if  execution  should  be  terminated.  Under  this  run  type,  if 
modifications  were  read  then  the  pressures  and  flows  are  again  calcu- 
lated and  printed.  Any  number  of  modified  data  sets  can  be  determined 
and  solved  in  this  manner. 
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2. 


RUN  1:  This  run  type  is  similar  to  the  one  above  except  that  no 

solutions  for  pressures  and  flows  are  produced  for  modified  data  read 
under  the  subroutine  INPUT2 . A solution  is  obtained  for  the  original 
data  read  under  subroutine  INPUT. 

3.  TEST:  In  this  run  type  data  are  read  and  tested  but  no  solutions  are 
produced.  Run  types  2 and  3 can  be  used  to  test  input  data  which  are 
to  be  used  later  under  another  run  type  and  under  a low  cost  priority. 

4.  FIND  F:  In  this  run  type,  the  rate  of  pressurization  air  necessary  to 

produce  a specific  pressure  differential  across  a shaft  is  determined. 

Figure  1 is  a flow  chart  for  the  main  program  logic  under  the  conditions 
of  a run  type  RUN. 

5 . 2 INPUT  Subroutine 

This  routine  reads  the  data  that  are  necessary  for  a flow  analysis  of 
the  stairwell  or  elevator,  including  an  analysis  of  the  rest  of  the  building. 
These  data  consist  of  the  following: 

1.  Outside  temperature. 

2.  Temperature  throughout  the  building. 

3.  Outside  pressures.  These  can  be  entered  or  calculated  as  described 
earlier. 

4.  Description  of  the  flow  network  including  flow  coefficients  and  flow 
exponents  for  all  connections  and  the  net  air  flows  to  each  space  due 
to  the  air  conditioning  system  or  due  to  a pressurization  system. 

The  data  above  can  be  entered  in  either  SI  units  or  in  engineering 
units.  There  are  two  modes  for  the  description  of  the  flow  network.  The 
first  is  a general  mode  which  allows  any  building  geometry  and  any  connec- 
tions between  building  spaces  that  the  user  chooses  to  define.  The  second 
description  mode  is  simplified  in  that  each  floor  is  represented  by  one 
building  space  and  the  flow  coefficients  and  exponents  are  all  the  same  in 
the  vertical  direction  unless  they  are  defined  by  exception  to  be  different. 
This  simplified  data  input  greatly  reduces  the  quantity  of  data  required  to 
define  a large  building.  Appendix  A contains  a detailed  description  of  both 
of  these  methods  of  input. 
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START 


INPUT  Subroutine 

Read  and  Print  Data 

Initialize  P Terms 
z 


CORR  Subroutine 
Calculate  Flow  Coefficients 
Corrected  for  Temperature 


Determine  if  Another 
Solution  is  Desired 
With  Modified  Data 


Figure  J . Flow  chart  for  main  program  logic  under  a run  type  RUN. 
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In  addition  to  reading  data,  this  subroutine  provides  temperature  and 
pressure  data  as  well  as  a complete  description  of  the  flow  network.  This 
routine  also  calculates  initial  estimates  of  the  hydrostatic  pressure 
differences.  When  data  is  entered  in  engineering  units  the  subroutine  UNITS 
is  called  which  converts  all  units  to  the  SI  system. 

5.3  CORR  Subroutine 

This  routine  calculates  flow  coefficients  corrected  for  the  temperature 
of  the  entering  air  using  eg.  (3.2).  Two  sets  of  these  corrected  coefficients 
are  calculated  for  each  flow  path  to  allow  for  flow  in  either  direction. 

5.4  INIT  Subroutine 

This  routine  calculates  initial  estimates  of  the  building  pressures  by 
a technique  used  by  Sander  [2] . In  this  technique,  mass  flows  are  considered 
linear  functions  of  differential  pressure  and  therefore  the  flow  equations 
can  be  expressed  and  solved  in  matrix  form.  In  this  estimate,  shaft  pressures 
are  considered  hydrostatic.  The  resulting  pressures  form  a starting  point 
for  the  iterative  solution  which  follows. 

5 . 5 BLDGP  Subroutine 

The  iterative  solution  for  the  building  pressures  and  flows  consists  of 
the  three  subroutines  BLDGP,  SHAFTP  and  PZAD.  The  subroutine  BLDGP  operates 
on  the  building  compartments  sequentially.  The  sum  of  all  the  mass  flows 
into  compartment  i is  calculated.  If  the  absolute  value  of  this  sum  is  less 
than  a convergence  limit  then  eq.  (3.3)  is  considered  satisfied  and  the 
computer  proceeds  to  the  next  compartment  or  returns  to  the  main  program. 
However,  if  the  absolute  value  of  the  sum  is  greater  than  the  convergence 
limit,  then  an  improved  estimate  of  the  pressure  at  compartment  i is  obtained 
by  the  regula  falsi  method  [16] . When  none  of  the  pressures  need  to  be 
modified  this  routine  passes  a convergence  signal  to  the  main  program. 

5 . 6  SHAFTP  Subroutine 

The  structure  of  this  routine  is  very  similar  to  that  of  BLDGP  except 
that  it  operates  on  shafts  sequentially.  The  sum  of  all  the  mass  flows  into 
shaft  i is  calculated.  If  the  absolute  value  of  this  sum  is  less  than  the 
convergence  limit  then  eq.  (3.4)  is  also  considered  satisfied  and  the  com- 
puter proceeds  to  the  next  shaft  or  returns  to  the  main  program.  However, 
if  the  absolute  value  of  the  sum  is  greater  than  the  convergence  limit,  then 
improved  estimates  of  the  shaft  pressure  are  calculated.  This  is  done  by 
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changing  the  pressures  at  the  bottom  of  the  shaft  and  then  recalculating  the 
shaft  pressure  by  eq.  (3.5)  . Again  the  regula  falsi  method  is  used,  and  if 
none  of  the  shaft  pressures  need  to  be  modified  a convergence  signal  is 
passed  to  the  main  program.  It  can  be  seen  from  figure  1 that  if  convergence 
is  achieved  in  both  BLDGP  and  SHAFTP,  then  the  subroutine  OUT  will  print  the 
solution.  Otherwise,  the  hydrostatic  pressure  differences  are  adjusted  in 
the  subroutine  PZAD. 


5.7  PZAD  Subroutine 

This  routine  calculates  hydrostatic  pressure  differences  by  eq.  (3.6) 
using  the  most  recent  pressure  estimates. 

5.8  OUT  Subroutine 

This  routine  outputs  mass  flows  and  pressures  for  the  flow  network  as 
well  as  the  differential  pressures  across  each  shaft.  If  the  data  input  was 
in  engineering  units  then  the  subroutine  UNITS  is  called  to  convert  variables 
to  the  engineering  system  before  output. 

5.9  INPUT2  Subroutine 

This  routine  reads  modifications  to  the  flow  network.  The  modifications 
can  consist  of  changes  in  the  net  flow  to  a stairwell  or  elevator  pressurization 
system  or  of  changes  to  flow  coefficients  between  either  between  building 
spaces  or  to  the  outside.  After  the  modifications  are  read,  a solution  to 
the  new  flow  network  is  obtained. 

6.  FUTURE  DIRECTION 

In  the  future  it  is  planned  to  analyze  pressurized  stairwell  systems 
with  pressurized  vestibules  as  has  already  been  done  for  systems  without 
vestibules  [10] . In  addition,  the  use  of  this  program  is  a part  of  a large 
project  to  study  the  feasibility  of  using  pressurized  elevator  shafts  for 
fire  evacuation  of  the  handicapped.  The  further  development  of  this  program 
so  that  it  can  easily  be  used  as  a design  tool  is  being  considered.  However, 
this  development  would  depend  upon  the  determination  that  a real  need  exists 
for  such  a design  tool. 
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DATA  INPUT  DESCRIPTION 


APPENDIX  A.  DATA  INPUT  DESCRIPTION 


Data  input  consists  of  the  following  elements: 

1.  Initial  data 

2.  Building  heights 

3.  Temperature  profiles 

4 . Outside  pressure  profiles 

5.  General  flow  network  input 

6.  Abbreviated  flow  network  input 

7.  FIND  F data 

8.  INPUT2  data 

Each  of  these  input  elements  is  described  in  detail  in  the  following 
sections.  Elements  1 through  4 are  always  required.  The  flow  network  can 
either  be  defined  by  a general  method  of  input  using  element  5 or  by  an 
abbreviated  method  using  element  6.  The  virtue  of  the  abbreviated  method  is 
that  less  data  are  required.  In  the  abbreviated  method  each  floor  is  repre- 
sented by  one  space  and  flow  coefficients  and  net  flows  are  all  the  same 
vertically.  However,  the  abbreviated  flow  network  input  element  also  allows 
a number  of  specific  exceptions  to  this  simple  model. 

Element  7,  FIND  F data,  is  used  for  a FIND  F run  type.  Element  8, 

INPUT2  data,  is  used  to  modify  the  existing  data  for  another  analysis.  In 
the  following  sections  the  input  required  for  each  of  the  eight  data  elements 
is  described  in  detail.  Each  block  or  group  of  blocks  below  represent  an  input 
card.  Unless  otherwise  stated  these  cards  are  unformatted,  that  is  the 
numbers  do  not  have  to  be  placed  in  specific  columns  of  the  card.  However, 
separate  pieces  of  numerical  data  must  be  separated  by  one  or  more  spaces. 
Examples  of  input  data  are  provided  in  Appendix  B. 

1 . Initial  data 

run  type^  either  RUN,  RUN  1,  TEST,  or  FIND  F (col.  1-6  must  start  in 
col.  1) 


project  title  (col.  1-72) 


A1 


outside  unit  indication  summary  output 

temperature  (°C,°F)  (1  for  SI,  2 for  Eng)  (0  for  none,  or  file  number)"^ 


2 . Building  heights 

N. , no.  of  building 
levels 


input  parameter 
(either  0 or  1) 


If  input  parameter  = 0,  then  heights  for  each  building  level  are  to  be 
individually  inputed  as  follows: 


(1) 


h (2)  h ( 3) 


h(i) 


(V 


o o o 


where  is  the  height  of  the  center  of  level  i above  the  ground  (m,  ft) 


If  input  parameter  = 1,  then  the  following  card  must  be  entered. 

h,^.  distance 

between  floors  (m,  ft) 


3 . Temperature  profiles 

no.  of  temperature  profiles 


^All  net  flows  are  inputed  at  standard  conditions  of  21°C  (70°F)  and  1 
atmosphere . 

2 

The  user  must  assign  this  file  before  program  execution. 
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For  each  temperature  profile  the  following  data  must  be  supplied. 


no.  of  temp, 
points 


level  temperature  level  temperature 
no.  (°C,  °F)  no.  (°C,  °F) 


level  temperature 
no.  (°C,  °F) 


O 


4 . Outside  pressure  profiles 


N 


po 

no.  of  outside 
pressure  profiles 


input  parameter 
(either  0 or  1) 


If  the  input  parameter  = 0,  each  outside  pressure  profile  is  entered  as 
follows : 


P P P 

o(l)  o(2)  o(3) 


Po(i) 


o(Nh) 


where  pQ(i)  i-s  the  outside  pressure  at  level  i. 


If  the  input  parameter  = 1,  the  outside  pressures  are  calculated  and  the 
following  data  are  required. 


V 

wind 

velocity 


h 

height  at  which 
velocity  is  measured 


n 

wind 

exponent 


pressure  coefficients  for  each  pressure  profile 


CW(1) 


CW(2) 


CW(Npo) 


o o o o 
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5 . General  flow  network  input 


5 . 1 General  building  data 


N 

no.  of  building  input  parameter 

compartments  (0  signals  non- abbreviated  input) 


For  each  compartment  the  following  data  must  be  supplied. 


N 

■eC 

no.  of  con- 
nections to 
other  spaces 


N 

j-O 

no.  of  con- 
nections to 
the  outside 


Ff 

net  flow 
( l/s , cfm) 


temperature 

profile 

number 


compartment 

level 


For  each  connection  between  this  compartment  and  other  spaces  in  the  building 
the  following  connection  data  are  required. 

number  of  CQ  3 X 

connected  flow  coefficient"3  flow  exponent 

space 


For  each  connection  to  the  outside  the  following  connection  data  are  required. 


outside  pressure 
profile  number 


C 

flow  coefficient 


X 

flow  exponent 


5.2  General  shaft  data 


no  of  shafts 


3 

All  flow  coefficients  are  entered  for  standard  conditions  of  21°C  (70°F)  and 

1 atmosphere.  The  units  for  all  flow  coefficients  are  Us-1 (Pa) or  cfm 

(”H20)*  1 *"X. 
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For  each  shaft  the  following  data  are  required. 


shaft  title  (col  1-18) 


C 

shaft  flow 
coefficient 


bottom 

level  of  shaft 


top  level 
of  shaft 


temperature 

profile 

number 


For  each  level  of  this  shaft  the  following  data  are  required. 


no 


. of  connections 
to  other  spaces 


no 


. of  connections 
to  the  outside 


net  mass 
flow  (£/ s,  cfm) 


For  each  connection  between  this  level  of  the  shaft  and  other  building 
spaces  the  following  data  are  required. 

no.  of  other  CQ  X 

spaces  flow  coefficient  flow  exponent 


For  each  connection  between  this  level  of  the  shaft  and  the  outside  the 
following  data  are  required. 


outside  pressure 
profile  number 


flow  coe¥ficient 


X 

flow  exponent 


A5 


6 . Abbreviated  flow  network  input 


6 . 1 Specific  building  data 
N 

no.  of  floors  input  parameter 

in  building  (1  signals  abbreviated  input) 


The  following  data  are  entered  only  once  because  it  is  the  same  for  all 
floors  of  the  building  except  at  locations  where  specific  exceptions  are 
defined . 


N 

■eO 

no.  of  con- 
nections to 
the  outside 


no.  of  the 
temperature 
profile 


flow  c8efficient 
through  floors 


X 

flow 

exponent 
through  floors 


no.  of 
exceptions 


For  each  connection  to  the  outside  the  following  data  is  required. 

outside  pressure  C X 

profile  number  flow  coefficient  flow  exponent 


For  each  exception  the  following  data  is  required. 


exception  type 
(1,  2 or  3) 


building  level 
for  the  exception 


The  next  card  depends  upon  the  exception  type.  For  exception  type  = 1,  a 
net  flow  is  defined  for  this  building  level. 

Ff 

net  flow  for 

building  level  (£/s,  cfm) 
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For  exception  type  = 2 , an  exception  to  an  outside  connection  is  defined. 


no.  of  outside  C X 

pressure  profile  flow  coefficient  flow  exponent 


For  exception  type  = 3,  a new  compartment  is  added  at  the  building  level. 

In  this  case  the  following  cards  are  needed  to  define  this  new  compartment. 


C 

flow  coefficient 
between  new  compartment 
and  rest  of  the  floor 


X 

flow 

exponent 


no 


. of  connections 
to  the  outside 


Each  of  the  outside  connections  is  defined  as  follows. 

outside  pressure  C X 

profile  number  flow  coefficient  flow 

exponent 


6 . 2 Specific  shaft  data 


no.  of  shafts 


For  each  shaft  the  following  data  are  required, 
shaft  title  (col  1-18) 


Cg  bottom  level  top  level  temperature 

shaft  flow  of  shaft  of  shaft  profile  number 

coefficient 
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The  following  data  are  entered  only  once  because  it  is  the  same  for  each 
level  of  the  shaft  except  at  locations  where  specific  exceptions  are  defined. 

C flow  no.  of 

net  mass  flow  coefficient  exponent  exceptions 

flow  to  building 

( l/s , cfm)  » 


For  each  outside  connection  the  following  data  are  required. 

outside  pressure  C X 

profile  number  flow  coefficient  flow 

exponent 


N 

rO 

no.  of  con- 
nections to 
the  outside 


For  each  exception  the  following  data  are  required. 


exception  type 
(1,  2 or  3) 


level  of 
shaft 


The  next  card  depends  upon  the  exception  type.  For  exception  type  =1,  an 
exception  is  defined  to  the  net  flow  into  the  level  of  the  shaft. 

net  mass^flow 
( H/s , cfm) 


For  exception  type  = 2 , an  exception  is  defined  to  the  outside  connection 
for  this  shaft. 

outside  pressure  C X 

profile  number  flow  coe¥ficient  flow 

exponent 
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For  exception  type  = 3 , an  exception  is  defined  to  the  flow  coefficient 
between  the  shaft  and  the  building. 

C 

flow  coefficient 


7 • FIND  F input. 

For  a FIND  F run  type  the  following  data  is  required. 


shaft  level 

number  of  shaft 


desired  differential 
pressure  across 
shaft 


convergence 

limit 


8.  INPUT2  data 


This  routine  is  used  to  modify  the  existing  data  for  a rerun. 


N 

no.  of 

modifications 


If  Nm  = 900,  program 
execution  is  stopped. 


For  each  modification  the  following  data  are  entered. 


modification 

type 


IS 

Integer 


B 

real 

number 


For  modification  type  = 1,  the  net  flow  (F^  = B)  at  each  level  of  shaft,  IS, 
is  changed. 

For  modification  type  = 2,  the  flow  coefficient  (C  = B)  between  each  level 
of  shaft  IS  and  the  building  is  modified. 

For  modification  type  = 3,  the  flow  coefficient  (CQ  = B)  to  the  outside  is 
changed  for  each  outside  connection  from  shaft  IS.  However,  if  IS  = 0,  then 
the  flow  coefficient  for  each  outside  connection  from  all  building  compart- 
ments is  modified. 
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For  modification  type  = 4,  the  flow  coefficient  (Cq  = B)  between  two  spaces 
(IS  and  JS)  is  changed.  In  addition  to  the  modification  type,  IS  and  B the 
second  space  number  is  entered  on  the  next  line. 


JS 

second  space 
number 
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APPENDIX  B.  INPUT  EXAMPLES 


APPENDIX  B.  INPUT  EXAMPLES 


1.  Example  1 

A ten  story  building  with  a pressurized  stairwell  and  no  vertical  leakage 
within  the  building  is  heated  to  70°F  when  the  outside  temperature  is  -20°F. 
The  stairwell  temperature  is  60°F  at  the  tenth  floor  and  50°F  at  the  bottom 
floor.  The  stairwell  is  pressurized  by  a net  400  cfm*  per  floor.  The  wind 
is  30  mph  at  a height  of  30  ft  and  the  wind  exponent  is  0.14.  This  building 
has  connections  to  the  outside  in  two  directions.  The  flow  coefficients  are 
the  same  vertically  and  are  listed  in  table  Bl  and  the  flow  exponents  are 
expected  to  be  0.5. 


Connection 

Cq/  flow  coefficients* 

location 

cfm  per  /"^O 

Between  stairwell  & 
building 

900 

Between  building  & 

outside  into  the  wind 

800 

Between  building  & 
outside  away  from 
the  wind 

800 

Table  Bl.  Flow  coefficients  for  example  1. 


1.1  General  Input 

In  the  general  flow  network  input  the  building  compartments  can  be  located  on 
any  floor  and  there  can  be  any  number  of  compartments  per  floor.  In  this 
case  there  is  only  one  compartment  per  floor  so  the  floor  levels  are  chosen 
to  agree  with  the  compartment  numbers. 

Therefore,  space  1 through  10  refer  to  floors  1 through  10  and  spaces  11 
through  20  refer  to  the  spaces  in  the  stairshaft.  Connections  between  two 
spaces  only  need  to  be  defined  for  one  of  the  spaces.  In  this  example  the 
connections  between  the  building  and  the  shaft  are  all  defined  for  the  shaft 
spaces . 


* 

At  standard  conditions 


Bl 


RUN 

Initial 

EXAMPLE  OF 

GENERAL 

FLOW 

NETWORK  INPUT 

data 

-20 

2 0 

building 

10 

1 

height* 

5 

10 

te.mpzAatu.sie 

2 

pAofille* 

1 

1 70 

2 

1 60 

1 

50 

out*lde 

2 

1 

pAe**uAe 

30 

30  .15 

pAofille* 

0. 

7 - 

0.7 

General 

building 

data 

10 

0 

Ut 

0 

2 

0 

1 

1 

{^loOA 

1 

800 

.5 

2 

800 

.5 

1 nd  1 

0 

2 

0 

1 

2 

1 \looA 

1 

800 

.5 

2 

800 

.5 

3 Ad  1 

0 

2 

0 

1 

3 

filoOA 

1 

800 

.5 

2 

800 

.5 

4th 

0 

2 

0 

1 

4 

filOOA 

1 

800 

.5 

2 

800 

.5 

5th 

[ 0 

2 

0 

1 

5 

filooA 

1 

800 

.5 

l 2 

800 

.5 

6th 

1 0 

2 

0 

1 

6 

t \l00A 

, 1 

800 

.5 

I 2 

800 

.5 

7th 

0 

2 

0 

1 

7 

{ \l00A 

1 

800 

.5 

2 

800 

.5 

Sth 

0 

2 

0 

1 

8 

t \l00A 

1 

800 

.5 

2 

800 

.5 

9th 

0 

2 

0 

1 

9 

faloOA  i 

1 

800 

.5 

2 

800 

.5 

1 Oth 

0 

2 

0 

1 

10 

ilOQA 

1 

800 

.5 

8 

800 

.5 

B2 


1 


I 

Gumnal 
ihCL^t 
data 

STAIRWELL 

| 72000  1 10  2 


u t 

1 

0 

400 

llooh. 

1 

900 

.5 

2 nd 

1 

0 

400 

hloon 

1 

900 

.5 

3 A.d 

1 

0 

400 

hloon  \ 

1 

900 

.5 

4th. 

1 

0 

400 

6 Iooa 

1 

900 

.5 

5th  1 

filoot i 1 

x 

0 

400 

1 

900 

.5 

6th 

1 

0 

400 

hloon 

1 

900 

.5 

7th 

1 

0 

400 

iloofi 

l1 

900 

.5 

8th 

1 

0 

400 

faloon. 

1 

900 

.5 

9th 

1 

0 

400 

filooh 

1 

900 

.5 

10  th 

1 

0 

400 

faloon 

1 

900 

.5 

znd 

data 

j 900 

1.2  Abbreviated  Input 


This  problem  can  also  be  solved  by  use  of  the  abbreviated  flow  network  input. 
This  reduces  the  quantity  of  data  required  and  to  that  extent  reduces  the 
chance  of  input  error.  The  abbreviated  data  is  as  follows: 


RUN 

Initial  EXAMPLE  OF  ABBREVIATED  FLOW  NETWORK  DATA 
data  -20  2 0 


building  ( 10  1 

height*  { 5 10 


te.mp<inatune. 
pn.o  kll<i.t> 


2 

1 1 70 

2 1 60  1 50 


B3 


outAtdz 

pHZt>bUHZ 

psiofitZeA 

ab  bazvtatzd 

bui.ldi.ng 

data 


abba2.ui.atzd 

Akafit 

data 


2 1 

30  30 

0.7  -0.7 

10  1 

I 2 1 

i 1 800 

1 2 800 


.15 


0 .5 

.5 

.5 


!1 

STAIRWELL 
72000  1 10 

0 400  900  .5 


end  { 900 

data 


0 


2 

0 


1.3  Example  1 Output 

The  output  for  this  example  is  the  same  for  both  these  forms  of  input.  This 
example  output  is  given  in  appendix  C. 


2 . EXAMPLE  2 


This  is  a 30  story  building  which  is  70°F  inside.  The  outside  the  air 
temperature  is  -5°F  and  there  is  no  wind.  This  building  has  a stairwell  and 
an  elevator.  The  flow  coefficients  which  are  generally  the  same  vertically 
are  listed  in  table  B2  and  the  flow  exponents  are  taken  to  be  0.5. 


Connection  location 

C flow  coefficient* 
o 

cfm  per  /"H ^0 

Between  stairwell  & 

900 

building 

Between  building  & 

1600 

outside 

Between  elevator  & 

1600 

building 

Table  B2 . Flow  coefficients  for  example  2. 


* 

At  standard  conditions 


B4 


In  this  example  a total  of  four  separate  cases  are  analyzed.  These  cases 
are: 

1.  First  floor  outside  door  open. 

2.  First  floor  outside  door  closed. 

3.  First  floor  door  between  the  elevator  and  the  building  open. 

4.  Stairwell  pressurized  at  a rate  of  900  cfm*  per  floor  and  all  doors 
closed. 

2.1  Example  2 - Input 

The  last  three  cases  are  handled  by  the  use  of  the  INPUT2  routine.  The  data 
are  as  follows: 


RUN 

tntttal  | 

EXAMPLE  INPUT  - 

- 30  STORY  BUILDING 

data 

-5 

2 

0 

butldtng  1 

30 

1 

kztghtA  1 

5 

10 

tnmp2.f1atan.11  j 

pn-ofitlH  1 

1 

1 

70 

outAtde. 

1 

1 

pn.HAAun.2. 

0 

1 

1 

data 

1 

abbn.Hvtatnd  | 

30 

1 

butldtng 

1 

1 

0 

0 1 

data 

1 

1600 

.5 

2 

1 

[ 1 

25000 

.5 

2 

Abbn.2.vtatnd  Ahafit  1 

STAIRWELL 

Shafit  Vata 

72000 

1 

30  1 

o 

0 

900 

.5  0 

Ahafit  2 

ELEVATOR 

80000 

1 

i — i 

o 

CO 

0 

0 

1600 

.5  0 

Cai>2.  2 : 

mo  dtfitcatto  n 

1 

type.  3 

3 

0 

1600 

* 

At  standard  conditions 
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Col6<l  3: 

modtfitcatton 

1 

typo,  4 

4 

1 

25000 

31 

Ca-ae  4: 

modtfitcatto  n 

2 

typ<U>  1 and  1 

2 

1 

900 

1 

1 

400 

end 

( 900 

data 

2 . 2 Example  2 Output 


The  output  for  example  2 case  1 (the  data  above  not  including  modifications 
for  Case  2,  3 and  4)  is  given  in  appendix  C. 
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APPENDIX  C.  EXAMPLE  OUTPUT 


APPENDIX  C.  EXAMPLE  OUTPUT 


I 


X 

£ 

eg 

in 

in  eg 

O' 

O 

N eg 

<0 

O' 

ffi  eg 

ffi 

eg  eg 

N ffi  eg 

O' 
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PROGRAM  LISTING 


APPENDIX  D.  PROGRAM  LISTING 

* MAIN  PROGRAM 


S>NBS*PL  IBS.SHOW 

C 


A. MAI N 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


COMPUTER  PROGRAM  FOR  AIR  FLOW  ANALYSIS  IN  BUILDINGS 
SPECIFICALLY  FOR  ANALYSIS  OF  PRESSURIZED  STAIRWELLS 
AND  PRESSURIZED  ELEVATORS 


PROGRAM  VARIABLES 

C FLOW  COEFFICIENT  BETWEEN  BUILDING  POINTS 

CO  FLOW  COEFFICIENT  TO  OUTSIDE 

CS  FLOW  COEFFICIENT  OF  SHAFT 

E LIMIT  WITHIN  WHICH  CONVERGENCE  IS  ACCEPTABLE 

F NET  FLOW  INTO  POINT  I 

FC  FLOW  BETWEEN  INTERNAL  POINTS 

FF  FIXED  FLOW  INTO  POINT  I 

FO  FLOW  TO  OUTSIDE 

FSS  NET  FLOW  INTO  SHAFT  IS 

H HEIGHT  FROM  GROUND  TO  MIDPOINT  OF  FLOOR 

I BUG  OUTPUT  VARIABLE 

ICONV  INTEGER  USED  IN  SUBROUTINES  BLDGP  AND  SHAFTP 

IF  ICONV  = 0 THEN  THE  PRESSURES  WERE  UNCHANGED 
I FLOOR  FLOOR  LEVEL  WHERE  POINT  IS  LOCATED 
IT  POINTER  TO  TEMP  PROFILE  FOR  POINT  I 

ITS  POINTER  TO  TEMPERATURE  PROFILE  OF  SHAFT 

JC  POINT  NO*  CONNECTED  TO  POINT  I 

JOC  DIRECTION  OF  OUTSIDE  CONNECTION 

N NO*  OF  BUILDING  COMPARTMENTS 

NC  NO.  OF  INTERNAL  POINTS  CONNECTED  TO  POINT  I 

NCO  NO.  OF  OUTSIDE  CONNECTIONS 

NFS1  BOTTOM  FLOOR  OF  SHAFT 

NFS2  TOP  FLOOR  OF  SHAFT 

NH  NO.  OF  FLOORS 

NPO  NO.  OF  OUTSIDE  PRESSURE  PROFILES 

NS  NO.  OF  SHAFTS 

NS1  I VALUE  FOR  START  OF  SHAFT 

NS2  I VALUE  FOR  END  OF  SHAFT 

NT  TOTAL  NO.  OF  POINTS  (BLDG  AND  SHAFT) 

NTP  NO.  OF  TEMPERATURE  PROFILES 

P PRESSURE  AT  POINT  I 

PFO  OUTSIDE  PRESSURE  PROFILES 

PO  OUTSIDE  PRESSURE 

PS  PRESSURE  PROFILE  OF  SHAFT  - WORKSPACE 

PZ  PRESSURE  DUE  TO  ELEVATION  DIFFERENCE 

T TEMPERATURE  PROFILE  ARRAY 

TITLE  PROJECT  TITLE 

TITSH  SHAFT  TITLE 

X FLOW  EXPONENT  TO  INTERIOF  POINTS 

XO  FLOW  EXPONENT  TO  OUTSIDE 


c 

PROGRA 

M PARAMETERS 

c 

MB 

MAX 

NO. 

OF 

BUILDING  COMPARTMENTS 

c 

MM 

MAX 

NO. 

OF 

POINTS 

c 

MS 

MAX 

NO. 

OF 

SHAFTS 

c 

MC 

MAX 

NO. 

OF 

CONNECTIONS 

FOR  ANY  POI 

NT 

c 

M PO 

MAX 

NO  • 

OF 

OUTSIDE  PRESSURE  PROFIL 

ES 

c 

M TP 

MAX 

NO. 

OF 

TEMPERATURE 

PROFILES 

Dl 
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MAIN  PROGRAM 


C MFL  MAX  NO*  OF  FLOORS 


PARAMETER  MM =2 20*  MS= 2 • MC=4 • MP0=2 . MTP=2 . MFL=1 05 , MB  = 1 05 
COMMON  NT « P( MM) « C(MM,MC)  ,NC( MM) , JC{ MM.MC) . ITS(MS ) * 

1 FC ( MM«  MC )cPZ!MM»MC)»POCMM,  MPO ) ,CO(  MM*MPO)  , F ( MM  ) , PFO!  MF L,  MPO)  «, 

2 FF(MM) «FO( MM, MPO)  ,CS( MS) ,PS<  MFL ) .NS1 ( MS)  . NS2(MS) , 

3 FSS(MS) e N*  NS  » NPO  » I CON V • E . I BUG , X ( MM , MC ) ,XO(MM,MPO),TITSH(MS.3)r 

4 NH,H<  MFL)  . I FLOOR! MM) , T ( MTP • MFL ) . NF SI ( MS)  * NFS2 (MS) . IT (MB) #NTP 

5 9 NCO(MM)e JOC ( MM* MPO) e TOUT 
DOUBLE  PRECISION  P»PO,PS 
COMMON  /RUN/I RUN 

DATA  MRUN/6HRUN  / 

DATA  MR UN 1/6 HR UNI  / 

DATA  NTEST /6HTEST  / 

DATA  NFINDF/6HFIND  F/ 


READ  ITEST  TO  DETERMINE  RUN  TYPE 
I TEST  = RUN  FOR  NORMAL  RUN 

ITEST  = RUN 1 FOR  A NORMAL  RUN  OF  1ST  DATA  AND 
A TEST  OF  ALL  OTHER  DATA 
ITEST  =TEST  FOR  A TEST  RUN  OF  ALL  DATA 

ITEST  = FIND  F TO  FIND  FIXED  FLOW  » FF « FOR  A SET  PRESSURE 
ACROSS  A FLOOR  OF  A SHAFT 

RE  AD ( 5o 804 ) ITEST 
NITER=500 
I R U N = 1 

CALL  INPUT  TO  READ  DATA 
CALL  INPUT 

SET  UP  FOR  FIND  F RUN  TYPE 
IS  = SHAFT  NO* 

IF  = FLOOR  NO. 

DPS  = SET  PRESSURE  ACROSS  SHAFT  IS  AT  FLOOR  IF 
EDP  = ALLOWABLE  CONVERGENCE  LIMIT 

IF(  ITEST  *NE.  NF I NDF ) GO  TO  10 
READ! 5* 700) IS. IF. DPS. EDP 
I DP  = NS 1 ( I S ) + I F— 1 
JDR  = JC(  I DP . 1 ) 

I 0=0 

N 1 = NS  1 ( IS) 

N 2 = NS2 ( IS) 

WRITE! 6e 807) IS. IF, DPS, EDP 
0 E = 0 * 2 

0 I CS=1 

IF!  ITEST  • EQ*  NTE ST ) GO  TO  41 

IF!  I R U N • GT • 1 .AND*  ITEST  • EQ  * MRUN11G0  TO  41 
TEMPERATUPE  CORRECTION 
CALL  CORR 

CALL  I NT  T TO  INITIALIZE  PRESSURE  ARRAY  , P 


D2 


no  non  o n o n •£•  r>  o o r>  non  * n o o w oon  non 


* MAIN  PROGRAM 

C 

22  CALL  INIT 

C 

C 

C DO  LOOP  TO  30  IS  ITERATIVE  SOLUTION  TO  PRESSURE  ARRAY 

C 

24  DC  30  ITER=1, NITER 

C 

C CALL  BLDGP  TO  SOLVE  FOR  BUILDING  PRESSURES 

C 

CALL  BLOGP 
I CB= ICGNV 

I F ( I CB  • EQ  • 0 .AND.  ICS  .EQ.  0 » GO  TO  40 

CALL  SHAFTP  TO  SOLVE  FOR  SHAFT  PRESSURES 

CALL  SHAFTP 
I CS  = I CONV 

I F ( ICB  * EQ » 0 .AND.  ICS  .EQ.  0)G0  TO  40 

CALL  PZAD  TO  CALCULATE  PZ  TERMS 

CALL  PZAD 
CONTINUE 

IF  ROUTINE  FAILS  TO  CONVERGE  IN  NITER 
ITERATIONS  PRINT  ERROR  MESSAGE 

W R I TE (6c  BOO ) 

CONTINUE 

WRITE(6«  80 1 ) I TER 

GC  TO  50  FOR  FIND  F RUN  TYPE 
IF ( ITEST  .EQ.  NFINDFSGO  TO  50 


CALL  OUT  TO  OUTPUT  SOLUTION 
2 CALL  OUT 

I F ( I TE  ST  .EQ.  NFINDF) WR ITE( 6*806) I S »FM 
1 IF<  ITEST  .EQ®  MRUN)WRITE( 6. 305 ) 

CALL  INPUT2  TO  CHECK  IF  ANOTHER  RUN  IS  REQUIRED 
WITH  MODIFIED  DATA 

CALL  INPUT2 

IF( ITEST  • NE«  NF INDF ) GO  TO  20 

ZERO  FF  FOR  SHAFT  IS  FOR  1ST  ITERATION  OF  A NEW  FIND  F RUN  TYPE 


00  44  I=N1,N2 

F F ( I ) = 0 * 

4 CONTINUE 

1 Q = 0 

GO  TO  20 


D3 


* 

c 

c 

50 


C 

c 

c 

52 


51 


53 

54 


56 


C 

C 

c 

58 

59 


60 


C 

c 

c 

c 

c 

c 

700 

800 

801 

804 

805 

806 
807 


808 

809 

900 


MAIN  PROGRAM 

SELECTION  OF  NEW  FIXED  FLOW  • FM  . FOR  FIND  F RUN  TYPE 

DPM  =P ( I DP )-P( JDP)  — DPS 
I Q- IQ+1 

I F ( 1 Q .GT.  30 ) GO  TO  42 
I F ( I BUG  • ECU  0 ) GO  TO  52 

WRITE!6»901 )D  PM  » DP  1 ,DP2 * F 1 ,F2 • I DP , J DP 

CHECK  CONVERGENCE  OF  DIFFERENTIAL  PRESSURE  , DPM 

I F ( AB S ( DPM ) • LT • EDPJGO  TO  42 
I F ( I Q .NE.  1 ) GO  TO  51 
DPI =DPM 
F1  = 0. 

F M = 500 • 

GO  TO  58 

IF( IQ  .NE.  2 ) GO  TO  53 
I F ( DPM  .LT.  0 • > GO  TO  42 
GO  TO  54 

IF(DPM*0P1  .GT.  0 • ) GO  TO  56 

DP2=DPM 

F2  = FM 

F M=  SORT  < (F1*F1  ) -DP  1 *( F 2 *F 2-F 1 *F 1 )/( DP2-DP1  ) ) 

GO  TO  58 
DPI =DPM 
FI  =FM 

FM=SQRT! ! F1*F1  )-DP 1 *!  F2*F  2-F  1 *F  1 )/( DP2-DP1  ) ) 

ASSIGN  FM  TO  FF  FOR  SHAFT  IS 

DO  59  I=N I . N2 
FF( I )=FM 

I F ( I BUG  .EQ.  0 ) GO  TO  60 
WRITE!6.900 )FM 
I CS=1 

IF( IQ  .GT.  3) GO  TO  24 
GO  TO  22 

END  OF  SELECTION  OF  MEW  FF  FOR  FIND  F RUN  TYPE 


FORMAT  STATEMENTS 
FORMAT!  ) 

FORMAT (/////5X, 35! 1H1 >//5X, 

+ 3 5H  F A I LURE  OF  MAIN  PROGRAM  TO  CONVERGE  //5X, 35 ! 1 HI  )// > 
FORMAT!  10X.I5.5X.il  HI TEP  AT  I O NS  ) 

FORMAT ! A6 ) 

FORMAT! 1 H 1 ) 

FORMAT! 10X.20HFIXED  FLOW  IN  SHAFT  .I2.3H  IS.F3.2/1H1) 

FOFMAT! 10! /) .10X.28HRUN  TO  DETERMINE  FIXED  FLOW  » 

1 18HPEP  FLOOP  IN  SHAFT  .I2/10X. 

2 41HSO  THAT  THE  PRESSURE  DIFFERENCE  AT  FLOOR  .13, 

3 3H  IS.F8.2.7H  PA SC AL/ 1 OX  , 2 8H W I TH  A CONVERGENCE  LIMIT  OF  # 

4 F8 • 5 o 7H  PASCAL# 10! /)  ) 

FORMAT! 10X# 10HNET  FLOWS  / /& ! 6X* 1HI .4X.3HNF  )/) 

F0PMAT!/10X»1 5HNET  FLOW  SHAFT  .I3.3H  IS.F6.1) 

FORMAT! 10X ,4HFM  =,F10.2) 
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901  F OR  MAT ( // 1 0X»  5HDPM  =* E 1 2. 4* 2X* 5H0P 1 =• E 1 2 • 4* 2X » 5HDP2  = *E12.4 

+ / 1 OX®  AHF 1 =9E 12. 4.2X* AHF 2 =. E l 2# 4 . 2X . 5HI DP  = « I 4 • 2X , 5H JOP  =.14) 
END 

®HDG9P  * SUBROUTINE  INPUT  .L . 1 
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* 


EDNB 

C 

c 

c 

c 


S*PLIB$.SHOW  A. INPUT 

SUBROUTINE  INPUT 

THIS  ROUTINE  READS  AND  PRINTS  DATA 
AND  INITIALIZES  PZ  ARRAY 

PARAMETER  MM=220 , MS=2,MC=4 ,MP0=2 * MTP=2 , M FL=1 05 . M B= 1 0 5 
COMMON  /PZZ/  PGZ 

COMMON  /I 0/T I T LE ( 12)  , I OUT, IUN I T , NCOMP ( MFL ) , S NCOMP ( MFL ) 
COMMON  NT  o P ( MM ) , C ( M M , MC ) , NC ( MM ) , JC ( MM , MC ) , I T S( MS ) , 

1 FC (MM* MC) e PZ ( MM, MC) ,PO(MM,MPO) ,CO( MM, MPO) , F ( MM ) , PFO( MFL i 

2 FF (MM) « FO( MM, MPO) ,CS( MS) ,PS( MFL ) ,NS1 ( MS) , NS2( MS) • 

3 FSS(MS) o N,NS ,NPO , I CONV , E , I BUG , X C MM , M C ) ,XQ(MM,MPO),TITSH« 

4 NH,H( MFL ) , IFLOOR  C MM) ,T( MTP , MFL ) , NF S 1 ( MS)  , NFS2( MS).  IT (MB] 

5 c NCO( MM) c JOCC MM, MPO) . TOUT 
DOUBLE  PRECISION  P,PO,PS 

DIMENSION  II (MFL) , TT ( MFL ) , P AR < 7 ) , C W ( MPO ) , P H( MFL ) ,NZZ(MM) 
DATA  PAR/3H  MM , 3H  MS,3H  MC ,3HMP0 , 3HMT P, 3HMFL , 3 H MB/ 

DATA  IBUG/0/ 

READ  AND  WRITE  PROJECT  TITLE 

READ ( 5P  600 ) (T  ITLE ( I ) , 1 = 1 , 12) 

W R I TE (6.601)  ( T I TLE ( I ) ,1=1  ,12) 


READ  GENERAL  DATA 


TOUT  = OUTSIDE  TEMPERATURE 
I UN I T = 1 FOR  SI  UNITS 
= 2 FOR  ENG  UNITS 
I OUT  = 0 FOR  NO  SUMMARY  OUTPUT 
OTHERWISE  I OUT  IS  FILE  NO.  TO 
WHICH  SUMMARY  OUTPUT  IS  WRITTEN 


READ( 5, 700 )TOUT , I UN  IT, I CUT 
WR I TE(6P4  1 1 )TOUT, I UNI T, I OUT 
I F ( IUN I T • GT • 2 .OR.  IUNIT  ,LT, 


1 )G0  TO  105 


READ  HEIGHTS 

NN=0  FOR  INPUT  OF  ALL  HEIGHTS 
NN=1  FOR  CALCULATION  OF  HEIGHTS 


89 


97 


RE 

AD  ( 5 

9 

700 

) 

NH,  NN 

WR 

I TE  ( 

6 

9 4 1 

2 

) NH  , 

NN 

I F 

( NH 

• 

LE. 

MFL  ) 

GO 

TO 

I P 

A P = 6 

GO 

TO 

1 

10 

I F 

( NN 

• 

EQ. 

1 ) GO 

T 

0 

97 

RE 

AD(  5 

s 

700 

) 

( H ( I 

) 

, 

I 

= 1. 

WR 

I T E ( 

6 

9 41 

3 

) ( H( 

I 

) 

e 

1=1 

GO 

TO 

9 

Q 

RE 

AC  ( 5 

0 

700 

) 

H ( 1 ) 

9 

D 

H 

WR 

I TE  ( 

6 

e 4 1 

4 

) H ( 1 

) 

, 

D 

H 

DO 

98 

I 

= 2. 

NH 

I M 

= 1-1 

i NH) 


MPO)  , 

MS, 3 ) „ 

« NTP 
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* SUBROUTINE  INPUT 

98  H( I )=H( IM)+DH 
C 

C R E AO  TEMPERATURE  PROFILES 

C 

99  READ! 5* 700  )NTP 
WRITE(6*415)NTP 

I F ( NTP  « LE  o MTP )G0  TO  90 
I PAR=5 
GO  TO  110 

90  DO  3 IP=1 e NTP 

READ! 5C  700 ) NNN  « ! I I ! J) ,TTI J) .J=l *NNN  J 
WRITE! 6o 416) NNN, ! I I < J ) * TT ( J ) * J = 1 ,NNN) 

IF( NNN  • GT • 1 ) GO  TO  2 
DO  1 IFF=1,NH 

1 T( IPc IFF )=TT I 1 ) 

GO  TO  3 

2 J = 1 
JP1  =2 

DO  4 I FF= 1 c NH 

T ( IP* IFF) =TT( J ) + ! TT! JP1 ) — TT ( J ) ) * ( IFF— I I ( J ) ) / ( 1 1 ( J P 1 )- I 1 ( J ) ) 

IF(  IFF  .NE.  1 1 ( JP1  ) ) GO  TO  4 

I F ( JP 1 .EQ • NNN ) GO  TO  4 

J=JP1 

JP1=J+1 

4 CONTINUE 

3 CONTINUE 

C 

c 

C READ  OUTSIDE  PRESSURE  PROFILES 

C NN=0  FOR  INPUT  OF  ALL  PRESSURES 

C NN=1  FOR  CALCULATION  BY  POWER  LAW 

C 

READ! 5o 700 ) NPO  » NN 

WRITE! 6c417)NPO»NN 

IF! NPO  *LE • MPO ) GO  TO  91 

IP4R=4 

GO  TO  110 

91  IF!NN  .EQ.  1 ) GO  TO  81 

C 

C READ  ALL  OUTSIDE  PRESSURES 

C 

DO  t 1 = 1*  NPO 

6 R EAD (5*  700 ) PGZ , ! PFO ( J , I ) , J=1 ,NH) 

WPITE!6o418)PGZ  , ! PFO! J , I ) , J=1 ,NH) 

GO  TO  85 
C 

C CALCULATE  OUTSIDE  PRESSURES 

C PATMOS  IS  ATMOSPHERIC  PRESSURE  !PA) 

C 

81  READ!5*700)VW»HW»XW,(CW!I  ) » I =1  ,NPO) 

WRITE!6*419)VW»HW»XW»!CW! I) *1=1 ,NPO ) 

IFIIUNIT  • EQ*  1 ) VW=VW*0.277e 
IFdUNIT  • EQ  • 2 )VW=VW*0.4470 

PATMOS=101 325. 

T OQ=TOUT  + 273. 

IFIIUNIT  .EQ.  2)T0C  = !T0UT  +4  60.  ) / 1 • 8 

PVA=176.4*VW*VW/T00 

Z=— 0 .03417/T00 
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I F ( I UNI T .EQ.  2)Z=0.3048*Z 
CWM=CW< 1 ) 

I F ( NPO  • EQ*  1)60  TO  212 
DO  211  1=1, NPO 

I F ( CW  < I ) .LT,  CWM) CWM=CWI I ) 

211  CONTINUE 

212  PGZ=PATMOS*EXP(HI NH ) *Z ) +C WM*PV A* ( I H C NH ) /H W ) ** ( 2 .*XW  ) )-100. 
DO  210  1=1 »NH 

PHI  I ) = P ATMOS* EXP { HI  I) *Z  ) 

210  CONTINUE 

DO  82  I = 1 e NPO 
DO  82  J=1 • NH 

PFOIJe  I )=PH(J)+CW(  I)*PVA*(IHIJ) /HW ) **  ( 2 • * X W ) ) -PGZ 
82  CONTINUE 

C 

c 

C EUILDING  DATA  INPUT 

C NFLS  = NO.  OF  FLOORS  IN  BUILDING 

C I F 1 = LOWER  FLOOR  IN  SERIES  OF  SIMILAR  FLOORS 

C IF2  = UPPER  FLOOR  IN  SERIES  OF  SIMILAR  FLOORS 

C NOC  = NO.  OF  COMPARTMENTS  PER  FLOOR 

C NZ  = NO.  OF  CONNECTIONS  TO  COMPARTMENTS  ON  SAME  FLOOR 

C NA  = NO.  OF  CONNECTIONS  TO  COMPARTMENTS  ON  FLOOR  ABOVE 

C 

85  1=0 

SNCOMPI 1 ) =0 • 

READI5.700) NFLS 
WR  ITEI 6»  420 )NFLS 
IFINFLS  .GT.  NH ) GO  TO  106 
7 READ! 5e 700) IF1 , IF2, NOC 

WRITE I6#400)IF1 , IF2.N0C 
I F ( IF 1 .GT.  I F2 ) GO  TO  107 
NCOMPI IF1 ) =NOC 
I FP= I F 1 +1 

SNCCMPI IFP)=SNCOMP( I FI  )+NOC 
DO  10  I Z= 1 c NOC 
1=1+1 

R E ADI  5c 700 ) N Z , N A , N NO , FF I I ) , IT! I ) 

WRITE  16*401) NZ. NA , NNO , FF ( I)  , IT I I ) 

NZZ I I ) = NZ 

NN=NZ+NA 

I FL  OOP ( I ) = IF 1 

IFINN  .LE.  MC)GO  TO  111 

IPAF=3 

GO  TO  110 

111  I F I NNO  .LE.  MPO  ) GO  TO  112 
I P A R = 4 

GO  TO  110 

112  IF(ITII)  .GT.  NTP  .OR.  ITII)  .LT.  1 ) GO  TO  102 
NCI  I )=NN 

IFINZ  .EQ.  0 ) GO  TO  63 

C 

C INPUT  CONNECTIONS  TO  COMPARTMENTS  ON  SAME  FLOOR 

C 

READ(5t 700) I JCI  I, J ) »C! I »J  ) »X( I , J ) » J = l ,NZ) 

WRITEI69402) 

WRITE(6®403)  (JC(I*J),CI  I, J)  ,XII  , J ) , J= 1 , NZ ) 

DO  62  J=1 g NZ 
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* SUBROUTINE  INPUT 

62  JC(  It J)=JC( It J )+SNCOMP(  IF  1 ) 

63  I F ( NA  .EQ.  0)GO  TO  8 
C 

C INPUT  CONNECTIONS  TO  COMPARTMENTS  ON  FLOOR  ABOVE 

C 

NP=NZ+1 

READ!  5*  700  ) ( JC<  I*  J ) • C ( I , J ) . X ( I , J ) . J =NP  . NA  ) 

W R I TE ( 6 o 4 04 ) 

KRITE( 6.  403 ) ( JC( I , J ) .C  < I , J)  .X ( I , J ) , J=NP,NA  ) 

00  66  J = NP « NA 

66  JC(  Is J )=JC ( It  J ) +NC  CMP ( I FI  )+SNCOMP( IF1 ) 

8 NCO  < I ) =NNO 

I F ( NNO  .EQ,  0 ) GO  TO  10 
C 

C INPUT  CONNECTION  TO  OUTSIDE 

C 

RE AD (5c  700)  (JOC(I,JJ),CO(I,JJ) ,X0< I , JJ ) . JJ  = 1 ,NNO> 
WRI TE (6,405) 

WRITE(6,403)( JOC( ItJJ)tCO(I.JJ).XO(ItJJ)eJJ=lcNNO) 
DO  9 J J = 1 o NNO 
J=J0C(I9 JJ) 

9 P0(  Ic JJ )=PF0(  IF1 , J ) 

10  CONTINUE 

I F ( I F 1 ,NE.  I F2 ) GO  TO  11 
I F { I F 1 .EQ.  NFLS) GC  TO  20 
GO  TO  19 
C 

C ASIGN  DATA  FOR  FLOORS  SIMILAR  TO  FLOOR  IF1 

C 

11  I FP  = I F 1 + 1 

DO  17  I FF= I FP  » I F2 

NCOMP(IFF)=NOC 
IFF  F=  I FF+ 1 

SNC  CMP( IFFP) =SNCOMP( IFF)  + NOC 
DO  16  I Z= 1 1 NOC 

1 = 1 + 1 

11=1  Z+SNCOMP ( I F 1 ) 

I FLOOR ( I ) = I FF 
FF ( I ) =FF ( I 1 ) 

I T ( I ) = I T(  II  ) 

NN=NC ( II ) 

NNO=NCO (II) 

NC ( I ) =NN 
NCO ( I ) =NNO 

I F ( IFF  ,NE.  NFL S ) GC  TO  23 
NN=NZZ(  II  ) 

NC ( I ) =NN 

23  I F ( NN  .EQ.  0 ) GO  TO  14 

DO  12  J = 1 » NN 
C(It J)=C( lie J) 

X ( I * J ) = X( II , J ) 

J C ( I « J ) = JC (II  * J)  + SNC0MP(IFF)-SNC0MP(IF1 ) 

12  CONTINUE 

14  I p( NNO  .EQ.  0 ) G O TO  16 

DO  15  J J= 1 » NNO 
JOC ( I s JJ) =JOC(  I 1 , J J ) 

J=JCC(  It  JJ ) 

C0(  It  J J )=CC( I 1 . JJ) 
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XO ( I o J J )=XO( I 1 » JJ ) 

15  PO<  Is JJ)=PFO(  IFF, J) 

16  CONTINUE 

17  CONTINUE 

18  I F(  IF2  .EQ.  NFLSJGC  TO  20 

19  CONTINUE 


GO  TO  7 
20  N=I 


N2  = N 

I F ( N • LE.  MB ) GO  TO  114 
I PAR=7 
GO  TO  110 


C 

c 

c 


SHAFT  DATA  INPUT 


114  R E AC  ( 5s 700 )NS 

I F ( NS  .LE.  MS ) GO  TO  113 
I PA  R=  2 
GO  TO  110 

113  DO  100  IS=1,NS 

RE AD (5*  603) (TITSH(  IS, I ) .1=1 .3 ) 

WRI TE(6e  406 ) ( T ITSH ( I S . I ) • 1 = 1 * 3 ) 

READ(5*700)CS(IS)  » NFS I ( IS)  .NFS2<  IS ) * ITS( IS  ) 
WRITE(6*4Q7)CS( IS) » NF  SI < IS) »NFS2( IS) • ITS< IS) 

N 1 =N2+ 1 

N2=N1+NFS2 (IS )— NFS1 ( IS) 

NS1 ( I S ) =N1 
NS2 ( IS) =N2 
I FF  = NF  S 1 ( IS)-1 

READ (5s  700)NN0,FFF, JCP.CC  »XX»NNN 
WRI TE( 6. 408)NN0,FFF, JCP.CC.XX .NNN 
I F ( NNO  .EQ.  0 ) GO  TO  21 

RFAD( 5*700) ( JOC ( N 1 » J ) »C0( Nl.J) .XO(NltJ) • J = 1 »NNO) 
WRITE(6*403)(JOC(N1*J)*CC(N1*J)  ,X0(  N1  » J ) . J=1 sNNO) 

21  DO  24  I =N 1 « N2 
NC< I )=1 

NCO ( I )=NNO 
I FF  = I FF+  1 
I FLOOR ( I ) = I FF 
J C ( I * 1 ) = JCP+SNCOMP ( IFF) 

C ( I * 1 )=CC 
X ( I * 1 )=XX 
FF ( I )=FFF 

I F ( NNO  .EQ.  0 ) GO  TO  24 
DO  22  J = 1 » NNO 
JJ=JOC(NlsJ) 

P0<  Is  J ) =PFO ( IFF* JJ) 

JOC ( I*  J)  = JJ 
C0(  Is J ) = C0(N1  , J ) 

22  X0(  Is  J ) =X0 ( N1 » J ) 

24  CONTINUE 


C 

c 

c 

c 

c 

c 

C 


EXCEPTIONS  TO  GENERAL  SHAFT  INPUT 
NNN  = NO.  OF  EXCEPTIONS 
KE  = 1 FOR  FF  EXCEPTION 
KE  = 2 FOR  OUTSIDE  CONNECTION 
KE  = 3 FOR  INTERNAL  CONNECTION 


DIO 


SUBROUTINE  INPUT 


I F ( NNN  ,EQ.  0 ) GO  TO  100 
DO  69  I K = 1 » NNN 
REAC(5«700)KE, IFF 
WRITE(6c409)KE, IFF 
I =N  S 1 ( IS) +IFF-NFS1 ( IS) 

I F ( KE  • EQ*  1 ) GO  TO  41 
I F ( KE  • EQ • 2 ) GO  TO  42 
I F ( KE  .EQ.  3 ) GO  TO  51 
GO  TO  104 

41  READ(5»700) FF ( I ) 
WRITE(6*410)FF( I ) 

GO  TO  69 

42  READ(5o700) J.CCC.XXO 
WR  ITE(6,405) 

WRITE (6, 403) J,CCO,XXO 
NNC=NCO( I ) 

I F ( NNC  .EQ,  0 ) GO  TO  44 
DO  43  K= 1 « NNC 

I F ( JOC ( It  K)  .EQ.  J ) GO  TO  46 

43  CONTINUE 

44  NJC=NNC+1 
NCO ( I )=NJO 

47  PO(  IeNJO)=PFO( IFF, J) 

JOC< I * N JO ) = J 
CO( I„NJO)=CCO 
XO ( I«  NJO)  = XXO 
GO  TO  69 
46  N JO  =K 

KK=K+ 1 

I F ( CCO  .NE.  O ) GO  TO  47 
N J0=NNC-1 
NCO ( I )=NJO 

I F ( N JO  .EQ.  0 ) GO  TO  69 
DO  49  K=KK,NNC 
KM=K-1 

P0(  Ic  KM )=P0( I, K ) 

JOC ( I * KM ) = JOC ( I , K ) 

C0(  I c KM ) =CO ( I ,K ) 

49  X0(  I«KM)=XO(  I ,K) 

GO  TO  69 

51  READ(5o700)JCP»CC»XX 
WR ITE(6»402) 

WRITE<6«  403) JCPeCC«XX 
J= JC°+SNCOMP( IFF) 

NN=NC( I ) 

I F ( NN  .EQ.  0 ) GO  TO  53 
DO  52  K= 1 , NNC 

I F ( JC ( T e K ) .EQ.  J)GO  TO  55 

52  CONTINUE 

I F( CC  .NE.  0.  ) G 0 TO  53 
WRITE (6, 520) IS.KE, IFF 
GO  TO  69 

53  N J=  NN+ 1 
NC ( I )=NJ 

54  JC(  I e N J ) = J 
C(  I t NJ )=CC 

X ( I «N J ) = XX  * 

GO  TO  69 
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55  NJ=K 

KK=K+1 

I F ( C C .NE.  O • ) GO  TO  54 
N J=  NN— 1 
NC ( I ) = N J 

I F ( N J • EQ.  0 ) GO  TO  69 
DO  61  K = KK  * NN 
KM=K-1 

JCC  Ic  KM)=JC( It K) 

C(  ISKM)=C( I ,K> 

61  X ( I « KM ) =X ( I , K ) 

69  CONTINUE 

100  CONTINUE 

NT=N2 

I F ( NT  .LE.  MM ) GO  TO  160 
I P A R= 1 
GO  TO  110 
C 

C PRINT  OUTSIDE  TEMPERATURE 

C 

160  WR ITE ( 6,  60 1 ) (TITLE!  I ) * 1 = 1 *12) 

IFCIUNIT  • EQ • 1 )WR ITE( 6,800 )TOUT 

IF(  IUNIT  .EQ,  2)WRITE( 6*5  00 ) TOUT 
IFCIUNIT  • EQ.  2 )TOUT=( TOUT-32. >/1.8 
TOUT=TOUT+273t 
C 

C PRINT  HEIGHT  AND  TEMPERATURE  PROFILES 

C 

I F ( I UN I T .EQ.  1 )WRITE( 6,811 )( IP, IP=1 ,NTP> 

I F ( IUN IT  .EQ.  2 ) WR ITE(6,511 ) ( IP, IP=1 ,NTP) 

WR  I TE (6,813) 

DO  30  I FF  = 1 , NH 

30  WRITE(6o312)H(I FF ) «(T( IP, IFF) , I P=1 , NTP ) 

C 

C CCNVERT  TEMPERATURES  TO  DEG  K 

C 

DO  33  I FF  = 1 , NH 
DO  33  I P= 1 , NTP 

IFCIUNIT  .EQ.  2)T( IP, IFF) =(T( IP, IFF)— 32. )/l.8 
33  T(  IP,  IFF)=T( IP, IFF ) + 273. 

C 

C PRINT  OUTSIDE  PRESSURE  PROFILES 

C 

IF(  IUNIT  .EQ.  1 ) GO  TO  79 
WRITE(6,514)(IP,IP=1,NP0) 

WRITEC6, 81 3) 

DO  76  IFF=1,NH 
DO  77  J = 1 » NPO 

77  PFO(IFF, J)=PFO( IFF,J)/248.8 

WRITE! 6, 51 5) H( IFF  ) , (PFO( IFF , J ) , J=1 , NPO) 

DO  78  J= 1 c NPO 

78  PFOC  IFF, J)=PFO(  IFF, J)*248.8 

76  CONTINUE 

GO  TO  83 

79  WRITE(6#814) ( IP, IP=1, NPO) 

WRITE(60813) 

DO  31  IFF=1,NH  , 

WRITE(6#815)H( IFF) ,(PFO(IFF,J) ,J=1,NP0) 
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31  CONTINUE 

C 

C CORRECT  FOR  CONNECTIONS  ONLY  INPUTED  ONCE 

C 

83  00  60  1=  1 « NT 

NN=NC ( I ) 

IFCNN  .EQ.  0 ) GO  TO  60 
DO  58  JJ=1.NN 
J=JC<  Ic  JJ ) 

IF( J .EQ.  0 ) GO  TO  58 
NNJ=NC ( J ) 

IFCNNJ  .EQ.  0 ) GO  TO  57 

DO  56  I A= 1 • NN J 

I F ( JC  ( J o I A ) .EQ.  I ) GO  TO  58 

56  CONTINUE 

57  NNJ=NNJ+1 

I F ( NN J .LE.  MCIGO  TO  59 
I P A R = 3 
GO  TO  110 

59  NC(J)=NNJ 
JC( Je NNJ )=  I 

C(  JcNNJ)=C( I*  JJ) 

X ( Jt NN J ) =X ( I 9 J J ) 

I F ( J .GT.  N .OR.  I .GT.  N)GO  TO  58 
PZ( JsNNJ) =-PZ<  I . J J ) 

58  CONTINUE 

60  CONTINUE 
€ 

C CORRECT  UNITS 

C 

I F ( I UN I T .EQ.  2 ) C A L L UNITS 
C 

C INITIALIZE  PZ  FOR  BUILD  COMPARTMENTS 

C 

87  DO  40  I = 1 9 N 

NN=NC ( I ) 

I F ( NN  .EQ.  0 ) GO  TO  40 
I A= IT ( I) 

I FI=IFLOOR( I ) 

DO  38  J J=1 »NN 
J = JC<  Ie  JJ) 

I F J = I FLOOR ( JO 
I F ( I F I .EQ.  I F J ) GO  TO  38 
I 8=  IT  ( J ) 

TEMPA=0.5*(T(IA,IFI)+T(IB.TFJ)) 

PZ<  I e JJ)=3A62.*<H(  IFJ )-H<  IF  I ) ) /TEMP A 
38  CONTINUE 

40  CONTINUE 

C 

C INITIALIZE  PZ  FOR  SHAFTS 

C 

DO  50  IS=1,NS 
N 1 = NS 1 { IS) 

N2=NS2 ( IS )-l 
I TT= I T S ( IS) 

DO  45  I =N 1 . N2 
I F I =1  FLOOR ( I ) 

I F J= IF  1 + 1 
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* 

TEMPA=0.5*<T( ITTe I F I ) +T ( I TT . I F J ) ) 

PZ<  I«  1 )=3462.*<  H(  IFJ )-H( IF I ) ) /TEMP A 
45  CONTINUE 

50  CONTINUE 

C 

C CHECK  SHAFT  CONNECTIONS 

C 

00  240  I S= 1 » NS 
N 1 =NS 1 ( IS) 

N2=NS2< IS) 

DO  239  I=N1»N2 
NN=NC( I ) 

I F ( NN  .EQ.  0 ) GO  TO  239 
DO  236  J = 1 » NN 
JJ  = JC<  It  J ) 

I F ( I FLOOR ( I ) • NE.  I FLOOPC J J ) ) GO  TO  103 

236  CONTINUE 

239  CONTINUE 

240  CONTINUE 
RETURN 

C 

C 

C DIAGNOSTIC  OUTPUT 

C 

102  WRITE (6.902) I , IT( I ) 

GO  TO  109 

103  W R I TE  ( 6c  90 3 ) 

GO  TO  109 

104  WP ITE ( 6«  904 ) 

GO  TO  109 

105  WR I TE ( 6e  905 ) 

GO  TO  109 

106  WR  ITE( 6*  90 6 ) 

GO  TO  109 

107  WRI TE ( 6e  907 ) 

GO  to  109 

110  WPITE(6c910)PAR( IPAR) 

C 

C PRINT  CORRECTED  BUILDING  DATA 

C 

109  WR  ITE (6e 940 ) 

DO  70  1=1. N 

NN=NC ( I ) 

I F ( NN  .GT.  0 ) GO  TO  180 

WRI TE ( 6o  94 1)1 , I FLOOR < I)eIT(I).FF(I) 

GO  TO  182 

180  WRITE! 60  942)1  .IFLOOP(I).IT(I).FF(I).JC(I.l).C(I.l).X(I.l) 
I F ( NN  .EQ.  1 ) GO  TO  182 

WRITE(60943)(JC(I.J)sC(I«J).X(I.J)»  J = 2.NN) 

182  NNC  = NCO ( I ) 

I F ( NNO  .EQ.  0 ) GO  TO  70 

WRI TEC  6,94  4)  ( JOC< I . J) ,C0( I , J)  ,X0( I . J) • J=1  .NNO) 

70  CONTINUE 

C 

C PRINT  CORRECTED  SHAFT  INPUT  DATA 

C 

DO  80  I S= 1 » N S 

WRITEC6.816) ( TIT  SH (IS«I)*I=1»3) 
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c 

400 

401 

402 

403 

404 

405 

406 

40  7 

408 

409 

410 

411 

412 

413 

414 

415 

416 

417 

418 

4 19 

420 

500 

511 

514 

515 

520 

600 

601 

603 

700 

800 
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WRITE! 6*  806) IS*  CS!  IS)*  ITS! IS) 

N 1 = NS 1 ( IS) 

N2=NS2 ( IS ) 

WRITE!6*807) 

DO  75  I =N 1 * N2 
NN=NC< I ) 

IF ( NN  • GT • 0 ) GO  TO  72 
WPITE!6»801 )I FLOOR ( I ) * FF ( I ) 

GO  TO  74 

WRITE! 6* 808) I FLOOR! I) • FF!  I ) , JC ( I • 1 ) • C ( I • 1 ) , X ( I • 1 ) 

IF! NN  .EQ.  1 )GO  TO  74 

WRITE!6»809)!JC!IeJ)*C!  I * J ) ♦ X ( I * J ) * J=2 1 NN) 

NNC=NCO! I ) 

IF! NNO  • EQ*  0 ) GO  TO  75 

WRITE!6»810)  (JOC! I , J ) , C 0! I » J ) » XC ! I «J) . J=1  .NNO) 

CGNTINUE 

CONTINUE 

STOP 

FORMAT  STATEMENTS 

FORMAT! 5X* 5HIF1  =,I3.7H,  IF2  =,I3,7H,  NOC  =.I3) 

F OR  MAT ! 5X»  4HN  Z =,I3,6H  NA  = . I 3 , 7H , NNO  = ,13, 6H.  FF  =,F8.1, 

+ 7H,  IT  = • I 3 ) 

FORMAT !5Xe  25HCONNECTION  ON  SAME  FLOOR  ) 

FORM  AT ! 5Xe  3HJ  =*I3*5H,  C =*F8.1*5H*  X =,F5.2) 

FOF MAT! 5X* 26HC0NNECT ION  TO  FLOOR  ABOVE  ) 

FORMAT !5X« 22HC0NNECT ICN  TO  OUTSIDE  ) 

FOP  MAT !5X*3A6) 

FORMAT !5X. 4HCS  =.F9.1,8H,  NFS1  =,I3,8H,  NFS2  =,I3,7H,  ITS  = ,I3> 
FORMAT !5X»  5HNN0  =,I3,7H,  FFF  =,F8.1*5H*  J = *I3*5H*  C = * F8,l* 

+ 5 He  X =*F5.2*7H*  NNN  = *13) 

FOPMAT (5X* 4HKE  =,13.  7H,  IFF  =.13) 

FORMAT! 5X« 4HFF  =,F8.1) 

FORMAT (5X* 6HT0UT  =,F6.0.9H,  IUNIT  =,I3.8H,  IOUT  =.I3) 

FOPMAT !5X« 4HNH  =*I3*6H«  NN  =*I3) 

FORMAT! 5X* 7HHEIGHTS  /!10F8#2)) 

FOPMAT  !5X* 6HH! 1 ) =,F8,2,6H,  DH  =,F8.2) 

FOPMAT !6X* 5HNTP  =.13) 

FORMAT!5X*  20HTEMPERATURE  PROFILE  /I  5 , ! 1 0 ! I 4 , F7 . 1 ) ) ) 

FORMAT  !5X,5HNP0  =,I3,6H,  NN  =,I3) 

F0RMAT!5X*5HPGZ  =*F12. 1/1 7HPRE  SSURE  PROFILE  /!10F12«1)) 

FOPMAT (5X« 4HVW  =*F6*1*6H*  HW  =.F6.1,6H.  XW  =.F4,2.6H.  CW  =. 

+ ! 1 0F4 • 2 ) ) 

FOPMAT !/5X. 6HNFLS  =,I3) 

FORMAT I//10X.20HOUTSIDE  TEMPERATURE  .F6.1.2H  F) 

FORMAT!  ///5X  « 6HHE I GHT , 5X • 2 9H TEMPERATURE  PROFILES  !DEG  F)  / 

+ 7X  ,2HFT, 3X*  191 6) 

F OP  MAT ! ////5X  * 6HHE I GHT  * 5X*  26HOUTSI DE  PRESSURE  PROFILES 
1 1 1 H ! IN  H20 ) /7X , 2HFT  .3X.8I10) 

FORMAT! FI 1.2.3X.8F10.3) 

FOF MAT!/// 5X . 1 5 HE PR OR  IN  SHAFT  . I 2 . 1 5HEXCE PT I ON  KE  = ,12. 

+ 2X.5HFL00R  ,13//) 

FORMAT! 1 2A6) 

FORMAT! 1H1///10X, 12A6///) 

FORMAT! 3A6* 12 ) 

FORMAT!  ) 

F OR  MAT ! //10X  « 20H0UTSIDE  TEMPERATURE  .F6.1.2H  C) 
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801  FORMAT! I13.F11.1 ) 

806  FORMAT!  1 OX  * 1 2H  SHAFT  NUMBER  . I 4/ 1 OX , 1 7H SHAFT  COEFFICIENT  .F10.1/ 

1  10X.20HTEMPERATURE  PROFILE  .14) 

807  F0RMAT(/21X.  5HF IX ED* 25 X# 4HFL0W. 12X.  4HFL0W/ 1 OX . 5HFL00R . 6X . 

1 4HFLQW.5X.1 2H CONNECTED  TO  • 6X . 1 1 HCOEFF I C IENT  . 6X . SHEXOONENT 

2 / ) 

8 08  FORMAT! I 13* FI  1.  1. 6X.  5HPOI NT . I 5* F 1 6. 1 • F 1 2. 2 ) 

809  FORMAT! 30  X.  5HPO INT  .I5.F16.1 .F12.2) 

810  FORMAT!30X.  7HOUTSI DE  . 13 • FI  6* 1 • F 12 • 2 ) 

811  FORMAT!  ///5X.6HHEIGHT.  5X  * 29HTEMPER ATURE  PROFILES  !05G  C)  / 

+ 7X.2HM  « 3 X. 1916) 

812  FOFMAT<F 1 1.2. 3X. 19F6. 1 ) 

813  FORMAT!/) 

814  F OR  MAT I////5X  » 6HHE I GHT  , 5X , 26HOUTS I DE  PRESSURE  PROFILES 

1 1 1 H < PASCALS ) /7X  « 2HM  .3X.8I10) 

815  FORMAT (FI  1 .2.  3X  » 8F 10*1) 

816  FORMAT! /// 10X  * 3A6  ) 

817  FOFMAT! 10X. 45HFL0W  COEFFICIENTS  CORRECTED  FOR  TEMPREATURE  ) 

902  FORMAT! 10! /)  , 10X. 1 1 HCOMPA RTMENT  .14/ 

1 10  X. 20 HTEMPER ATURE  PROFILE  .I4.17H  DOES  NOT  EXIST  / 

+ 10X* 16HPR0GRAM  STOPPED  ,10!/)) 

903  FORMAT! 10!/), 5X. 23HSHAFT  CONNECTION  ERROR  . 

1 /I  OX. 16HPR0GRAM  STOPPED  .10!/)) 

904  FORMAT! 10!/). 10X, 40HINPUT  ERROR  IN  EXCEPTIONS  TO  SHAFT  DATA 
1 /I  OX. 16HPR0GRAM  STOPPED  .10!/)) 

905  FORMAT! 10(/)o 10X. 37HINPUT  ERROR  IN  UNIT  TYPE  DESIGNATION  / 

1 1 OX. 16HPR0GRAM  STOPPED  .10!/)) 

906  FORMAT! 10!/). 10X. 37HINPUT  ERROR  NO.  OF  FLOORS  EXCEEDS  NH  / 

1 10X. 16HPR0GRAM  STOPPED  ,10!/)) 

907  FORM  AT (10(/)»10X*  2 5H INPUT  ERROR  IF1  .GT.  IF2  / 

1 10X. 16HPR0GRAM  STOPPED  ,10!/)) 

910  FOF  MAT! 10! /)» 10X.36HI NPUT  EXCEEDS  DIMENSION  PARAMETER  . A 3/ 

+ 10X. 16HPR0GRAM  STOPPED  .10!/)) 

930  FORMAT! 10X. 3A6) 

935  FORMAT!//  10X.26HFL0W  COEFFICIENTS  AS  READ  ) 

940  FOF MAT! 10X. 15HBUILDING  DATA  //34X * 1 1 HTEMPER ATURE  . 4X . 5HFIXED. 

1 1 2X, 2! 1 1 X,4HFLOW)/10X, 1 1 HCCMPA RTMENT  . 4X® 5HFL00R. 6X. 7HPR0F ILE. 

2 6X.4HFLOW. 5X. 1 3HCONNEC T I ON  TO  . 4X.  1 1 HCOEFF IC IE  NT  .4X, 

3 8HEXP0NENT  ) 


» I3.F1 1. 1.F13. 1 ) 


941 

FORMAT 

942 

FORMAT 

943 

FORMAT 

9 44 

FORMAT 

END 
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SUBROUTINE  CORR 


* 


3NBS*°L IBJ.SHOW  A. CORR 

SUBROUTINE  CORR 
C 

C THIS  ROUTINE  CALCULATES  FLOW  COEFFICIENTS 

C WHICH  ARE  CORRECTED  FOR  TEMPERATURE 

C 

°AR AMETEP  MM =2 20.  MS= 2* MC=4 , MP0=2 * MTP  = 2 • MFL= I 05 . M B= 1 05 
COMMON  /CORR/C1  ( MM. MC ) , C2 (MM ,MC  ) .CO  1 < MM.MPO) .C02( MM.MPO) 

COMMON  NT.  PC  MM)  ,C( MM.MC)  ,NC( MM) , JC( MM.MC) . IT3( MS  ). 

1 FCCMM, MC) .PZCMM, MC) , POCMM. MPO)  .CO( MM.MPO)  , F ( MM ) , PFOC  MF L, MPO)  . 

2 FF ( MM) .FO( MM, MPO)  ,CS(  MS)  ,PS < MFL ) • NS1 ( MS)  , NS2 ( MS) , 

3 FSS( MS ) . N. NS . NPO , I CON V. E , I BUG , X ( MM , M C ) , X O ( MM , MPO ) • TI T SH<  MS • 3 ) • 

4 NH.HC MFL ) , IFLOORC  MM)  ,T<  MTP.MFL) ,NFS1  ( MS)  . NFS2C  MS).  IT C MB) ,NTP 

5 . NCOC MM) . JOC( MM. MPO) .TOUT 

DOUBLE  PRECISION  P.PO.PS 
DO  12  1=1 . NT 

C 

CORRECT  C 
NN=NC ( I ) 

I F ( I .GT.  N ) GO  TO  1 
IP=  IT ( I ) 

GO  TO  4 
DO  2 I S=1 . NS 

I F ( I ,LE.  NS2 (IS)  .AND.  I • GE . NSl(IS))GO  TO  3 
CONTINUE 
WP ITE (6. 700 ) 

STOP 

IP=  IT  S ( IS) 

I FF  = I FLOOR ( I ) 

T 1 =T ( IP, IFF) 

I F ( NN  .EQ.  0 ) GO  TO  10 
DO  9 J=1 , NN 
J J = JC ( I. J ) 

Cl(I.J)=C(l9J)*((294./Tl)**<2.5*X(I.J)-0.75>) 

I F ( J J .GT.  N ) GO  TO  5 
I P=  IT ( J J ) 

GO  TO  8 
DO  6 I S = 1 9 NS 

I F ( J J .LE.  NS2 (IS)  .AND.  JJ  .GE.  NS1(IS))G0  TO  7 
CONTINUE 
WR  I TE ( 69  700 ) 

STOP 

IP= IT  S ( IS) 

IFF=IFLOOR( JJ) 

T2=T( IP, IFF) 

C2(I.J)=C(I.J)*((294./T2)**(2.5*X(I,J)-0.75)) 

CONTINUE 

CORRECT  CO 
C 

10  NNC=N  CO ( I ) 

I F ( NNC  .EQ.  0 ) GO  TO  12 
DO  11  J = 1 , NNC 

COl (I.J)=CO(I.J)*( (294./T1) **(2.5*XO(I.J)-0.75)) 

C02 (I  9 J ) = CO( I , J )*( ( 294 ./TOUT) ** (2.5*X0(  I , J )— 0.75) ) 

11  CONTINUE 
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* 


12  CONTINUE 
RETURN 

700  F0RMATC///10X, 36HPR0GRAM  STOPPED  IN  SUBROUTINE  CORR 
END 

3)  HOG  9 P * SUBROUTINE  INIT  .L. 


//) 
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SUBROUTINE  INIT 


* 


®NBS*PL IBS. SHOW  A. INIT 

SUBROUTINE  INIT 
C 

c 

C THIS  ROUTINE  INITIALIZES  THE  PRESSURE  ARRAY 

C 

PARAMETER  M M =2 20 . MS= 2 * MC =4  * MP0=2 . MTP=2 • MFL=1 05.MB=105 
PARAMETER  M8P=MB+1 

COMMON  NT.  P(MM).C(MM.MC).NC(MM). JC(MM.MC). ITS(MS). 

1 FC(MM.MC).PZ(MM.MC ) oPO(MM. MPO) *C0( MM. MPO) . F ( MM ) , PF 0{ MFL . MPO) . 

2 FF (MM) «F0( MM, MPO)  ,CS<  MS)  . PS ( MFL) .NS l ( MS ) . NS2(MS) . 

3 FSS(MS).N.NS , NPO . I CON V. E , I BUG . X ( MM , MC ) , X 0 ( MM , MPO)  . TI T SH ( MS . 3 ) . 

4 NH  » H ( MFL ) . I FLOOR ( MM) . T( MTP.MFL) .NFS1  ( MS)  . NFS2( MS) . IT (MB) . NTP 

5 . NCO< MM) . JOC ( MM. MPC) . TOUT 
DOUBLE  PRECISION  P.PO.PS 
DIMENSION  SC( MS) .SCO(MS) 

COMMON  /MAT/A(MB*M  BP ) . XX ( MB ) » NNN 
DOUBLE  PRECISION  A, XX 

NNN  = N 

e 

C CALCULATE  AVERAGE  OUTSIDE  PRESSURE 

e 

SUM  =0 • 

DO  10  J= 1 . NPO 
DO  10  1=1. NH 
10  SUM=SUM+PFO( I , J ) 

PA=SUM/(NPO*NH) 

C 

c 

c 

C THE  DO  LOOP  TO  STATEMENT  30  ESTIMATES 

C SHAFT  PRESSURES 

C 

DO  30  I S= 1 » NS 
C 

C CALCULATE  SHAFT  PRESSURE  DIFFERENCE  . DP 

C 

SUM  =0 • 

SUMN=0 • 

SUM  X=0 • 

NX  = 0 

N 1 = NS  1 ( IS) 

N2=NS2 (IS) 

DO  18  I=N 1 . N2 
SUM=SUM+FF< I ) 

NN=NC<  I ) 

I F ( NN  .EQ.  0.  ) GO  TO  16 
DO  15  J= 1 . NN 
SUMX=SUMX+X ( I . J ) 

NX=NX+1 

SUMN=SUMN+C( I. J) 

15  CONTINUE 
SC(  IS ) =SUMN 

16  NNC=NCO(I) 

I F ( NNO  .EQ.  0 ) GO  TO  18 
DO  17  J= 1 . NNO 
S UM  X = SUMX  + XO ( I . J) 
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N X=  NX+1 

SUMN=SUMN+CO< I « J ) 

17  CONTINUE 

SCO ( IS) =SUMN—  SC ( IS  ) 

18  CONTINUE 

A X = NX/SUMX 
DP2=SUM/SUMN 
S I GN= 1 • 

IP( DP  2 .LT.  0 • ) S I GN=— 1 • 

DP  = S IGN*(DP2**AX) 

C 

C CALCULATE  AVERAGE  TEMP  OF  SHAFT 

C 

SUM  =0 • 

IP=  ITS ( IS) 

DO  20  I=N1,N2 
IFF=IFL00R< I ) 

20  SUM=SUM+T ( IP, IFF) 

TA=SUM/(N2— Nl+1 ) 

C 

C ESTIMATE  PRESSURE  AT  BOTTOM  OF  SHAFT  • 

C 

HH=0.5*(H(NH)-H<1 ) )+H( 1 ) 

NF 1 =NF  S 1 < IS) 

P80T=PA+DP+3462.* ( HH-HINFl ) )/TA 

ESTIMATE  OTHER  SHAFT  PRESSURES 

P ( N 1 ) =PBOT 
NM=N2— 1 
DO  24  I=N1»NM 
IP1=I+1 

4 P ( I PI )=»< I )-PZ(  1,1) 

0 CONTINUE 

END  OF  SHAFT  PRESSURE  ESTIMATES 
C 

C SET  UP  MATRIX  FOR  BUILDING  COMPARTMENTS 

C 

NP 1 =N+ 1 
DO  50  1 = 1 o N 

NN  = NC(  I ) 

SUM  I 1 = 0. 

SUM  NP  = 0 • 

I F ( NN  .EQ.  0 • ) GO  TO  42 
DO  40  JJ  = 1 , N N 
J = JC(  I, JJ ) 

I F ( J .GT.  N ) GO  TO  34 
A ( I « J )=C(  I • J J ) 

SUM  I I = SUM I I-C(  I , JJ  ) 

SUMNP=SUMNP-C(  I , J J )*PZ ( I , J J ) 

GO  TO  40 

34  SUM  II=SUMI  I-C<  I , JJ  ) 

S UMNP  = SUMNP— C ( I » J J ) *P<  J) 

40  CONTINUE 

42  NNC=NCO( I ) 

I F ( NNO  .EQ.  0 ) GO  TO  46 
DO  45  K=1 , NNO 


PBOT 
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SUM  1 1 = SUM I I — CO ( I.K) 

45  SUMNP  = SUM NP—CO ( I tK)*PO(  I «K) 

46  A ( I « I ) = SUM I I 

A ( I . NP 1 )=SUMNP-FF(  I) 

50  CONTINUE 

C 

C WRITE  MATRIX 

C 

I F ( I BU  G • EQ • 0 ) GO  TO  84 
WR I TE ( 6 « 802 ) 

DO  52  1=1. N 

52  WRI TEC  6. 803) ( A< I . J ) . J=1 ,NP1 ) 

C 

C 

C CALL  ROUTINE  TO  SOLVE  FOR  INITIAL  BUILDING  PRESSURES 

C 

84  CALL  SIMEQ 
C 

C OUTPUT  INITIAL  PRESSURES 

C 

I F ( IBUG  .EQ.  0 ) GO  TO  89 
WRITE(6»  800 ) 

WRITE (6*801)(I. XX (I)*I=1.N) 

NN=NS 1 ( 1 ) 

WPITE(6.801 ) ( I»P( I ). I=NN»  NT) 

C 

C 

C ASSIGN  BUILDING  PRESSURES 

C 


89 

DO  90  1=1. N 

90 

P<  I ) =XX { I ) 
RETURN 

800 

FORMAT! ///8<6X • 1HI , 4X.3HP 

)/) 

801 

FORMAT (8(17. F7.1 ) ) 

802 

FORMAT ( ///l OX . 20HMATR IX 

COEFFICIENTS  /) 

803 

FOF  MAT ( 10X.  1 1 FI  1 . 1 ) 
END 

SHOGs 

P * 
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SUBROUTINE  BLDGP 


* 


$NBS*PL I B$. SHOW  A. BLDGP 

SUBROUTINE  BLDGP 
C 

c 

C THIS  ROUTINE  CALCULATES  STEADY  STATE  PRESSURES 

C FOR  BUILDING  COMPARTMENTS 

C 

c 

c 

PARAMETER  MM=220 * MS=2 * MC=4 , MPO =2 « MTP= 2 • MFL=1 05 . M B= 1 0 5 
COMMON  NT  p P( MM) #C(KM« MC) .NC(MM) . JC( MM.MC) , ITSCMS ), 

1 FC <MM« MC) « PZ (MM, MC) ,PO( MM, MPO)  ,CO<  MM, MPO)  , FC  MM ) , PF 0( MFL, MPO)  , 

2 FF ( MM) c FO( MM, MPO ) ,CS<  MS)  ,PS < MFL ) , NS 1 < MS ) ,NS2(MS), 

3 FSS ( MS ) * N, NS ,NPO, ICONV.E, I BUG , X ( MM , MC ) , X 0( MM , MPO ) , T I T SH ( MS , 3 ) , 

4 NH , H ( MFL ) ,IFLOOR( VM) , T< MTP , MFL ) ,NF SI (MS) ,NFS2(MS>* IT(MB).NTP 

5 c NCO( MM), JOC( MM, MPO) » TOUT 
DOUELE  PRECISION  P,PO,PS,PI 
I F ( I BUG  • GT • 0 )WRITE( 6,806) 

ITM=20 
I CONV=0 
DO  15  1=1, N 

CALCULATE  NET  FLOW  ,FI,  INTO  POINT  I 
F I = PFLOW ( I , P ( I ) ) 

CHECK  MAGNITUDE  OF  FI 
IF(AeS(FI)  .LT.  E ) GO  TO  15 
I C0NV=IC0NV+1 

SET  UP  PARAMETERS  FOR  ITERATION 
3D=1.0 
I PH  A SE  = 1 
DP  1=0. 

EE=0.2*ABS<  FI ) 

I F ( EE  .LT.  E ) EE=E 
S IGN=1 

I F ( F I .LT.  0.)SIGN=-1 

I K =0 

I F ( I BUG  .GT.  0) WRITE(6. 802) 

ITERATION  TO  RE  DU I CE  MAGNITUDE  OF  FN 
I K=  IK+1 

NEW  ESTIMATE  OF  PRESSURE  ,PI,  AT  POINT  I 
P I=P<  I )+SIGN*DP 

CALCULATE  NET  FLOW  , FN , INTO  POINT  I USING  PI 
FN  = FFLOW ( I , P I ) 

I F ( I BUG. GT.O) WRITE (6, 80 4) I , I K , F I , FN ,FP , DPI , DP , DPP , P I , I PH A SE 

C 

C CHECK  MAGNITUDE  OF  FN 

I F ( ABS ( FN  ) .LT.  EE ) GO  TO  10 
C 

C CHECK  NUMBER  OF  ITERATIONS 

I F ( I K .GT.  I TM ) GO  TO  25 

C 

C CHFCK  PHASE 
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IF( IPHASE  .EQ.  2 ) GO  TO  6 

C 

C CHECK  FOR  TRANSITION  FROM  PHASE  1 TO  PHASE  2 

I F( F I *FN  • LT • 0 • ) GO  TO  4 

ec 

C PHASE  1 

DPI=DP 
DP=5.0*DP 
F I = FN 
GO  TO  2 
C 

C °HASE  2 

4 IPHASE=2 

GO  TO  9 

6 IF(FI*FN  #GT.  O.IGO  TO  8 

C 

C NEW  DP  BETWEEN  DPI  AND  DP 

9 DPP=DP 
FP=FN 

DP  = DPI  4-  (DPP— DP  I >*FI/IFI-FN) 

GO  TO  2 
Cr 

C NEW  DP  BETWEEN  DP  AND  DPP 

8 FI=FN 

DPI =DP 

DP=DPI+ (DPP-DPI  )*FN/(FN— F°) 

GO  TO  2 

10  P(I)=PI 

15  CONTINUE 

C 

RETURN 

25  WPITE<6*800> 

STOP 

C 

C FORMAT  STATEMENTS 

C 

800  F OF  M AT ( /// 1 OX  * 20( 1 H* ) //✓! OX • 22HEXCESSI  VE  ITERATIONS  / 

+ 1 0 X*  8H IN  BLDGP  /// 10X * 20  < 1 H* ) //////) 

002  cOP MAT ( //I IX*  1HI*2X»2HIT*12X  »2HF I * 1 3X  « 2HFN  * 1 3X  * 2HFP  « 1 2X  * 3 HOP  I • 

■*•1  3X  • 2HDP  • 1 2X  * 3HDPP  « I3X.2HPI  «3X«5HPHASE  /) 

604  FOFMAT( 9X*  214, 3E1 5,4,4F15,6* 15) 

806  FORMAT!  /// 1 0 X, 6HBLDGP  ) 

END 

5)HDGo  R * SUBROUTINE  SHAFTP  ,L,1 
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SUBROUTINE  SHAFTP 
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<DNBS*°L  IB$,  SHOW  A.  SHAFTP 

SUBROUTINE  SHAFTP 

C 

C 

C THIS  ROUTINE  CALCULATES  STEADY  STATE  PRESSURES 

C FOR  SHAFTS 

C 

c 

PARAMETER  MM =220 , MS=2 , MC =4 , MPO=2 , M TP=2 , MFL =1 05  ,M8=105 
COMMON  NT,  P< MM)  ,C< MM ,MC)  ,NC(MM), JC< MM, MC)#  ITS(MS  )• 

1 FC(MM»MC) *PZ(MM*MC) • PO ( MM*  MPO ) , CO ( MM • MPO ),F( MM ) , PFO<  MFL,  MPO)  , 

2 FP (MM),FO(MM,MPO) ,CS( MS) . PS C MFL ) , N SI ( MS) . NS2(MS) . 

3 FSS( MS) • N*  NS ,NPO, ICONV,E.I8UG,X<MM,MC),XO(MM,MPO),TITSH<MS.3), 

4 NH,H<  MFL) * I FLO OR (MM) »T(MTP , MFL ) ,NF  SI ( MS)  , NFS2C MS  ) , IT  < M B ) , N TP 

5 ,NCO(MM), JOC(MM,MPO)*TOUT 
DOUBLE  PRECISION  P,PO,PS,PI 
I F ( I BUG  • GT  9 0 )WRITE(6,806> 

I TM  =2  0 
I C0NV=0 
DO  15  1=1, NS 

CALCULATE  NET  FLOW  ,FI,  INTO  POINT  I 
N 1 = NS  1 ( I ) 

FT  = SFLOW( I » P ( N 1 ) ) 

CHECK  MAGNITUDE  OF  FI 
IF( ABS(FI ) ,LT.  E ) GO  TO  15 
ICCNV=ICONV+l 

SET  UP  PARAMETERS  FOR  ITERATION 
DP= 1 • 0 
I PH ASE  = 1 
D P I =0 • 

EE=0.2*ABS(FI ) 

I F( EE  .LT,  E)EE=E 
S I GN= 1 

I F ( F I , LT • 0.)SIGN=— 1 
IK=0 

I F ( I FUG  • GT • 0) WRITE(6,802) 

ITERATION  TO  REDUICE  MAGNITUDE  OF  FN 
I K=  IK+1 

NEW  FST 1 MATE  OF  PRESSURE  ,PI,  AT  BOTTOM  OF  SHAFT  I 
P I=P(N1  )+SIGN*DP 

CALCULATE  NET  FLOW  , FN , INTO  SHAFT  I USING  PI 
FN=SFLOW< I , P I ) 

IF<  IBUG. GT.O ) WRITE < 6,804)  I,  I K , F I , FN , FP , DPI  , DP , DPP , P I , I PHASE 

CHECK  MAGNITUDE  OF  FN 
I F ( ABS ( FN ) *LT • EE ) GO  TO  10 

CHECK  NUMBER  OF  ITERATIONS 
IFC1K  , GT,  I TM ) GO  TO  25 

CHECK  PHASE 
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I F ( IPHASE  .EQ.  2 ) GO  TO  6 

CHECK  FOR  TRANSITION  FROM  PHASE  1 TO  PHASE  2 
I F( F I *FN  ,LT*  0 • ) GO  TO  A 

PHASE  1 
DPI=DP 
0P=5.0*DP 
F I =FN 
GO  TO  2 
C 

C PHASE  2 

4 I PH  A SE  = 2 

GO  TO  9 

6 IF<FI*FN  .GT.  0 . ) GC  TO  8 

C 

C NEW  OP  BETWEEN  DPI  AND  DP 

9 DPP=DP 
FP  = FN 

DP=DPI+(DPP-DPI )*FI/(FI-FN) 

GO  TO  2 
C 

C NEW  DP  BETWEEN  DP  AND  DPP 

8 F I = FN 

□PI =DP 

D P= DPI  + ( DPP— DP I )*FN/(FN-FP) 

GO  TO  2 

10  N2=NS2 ( I ) 

DO  11  I F =N  1 1>  N2 
I 1= IF+1-N1 

11  P ( I F ) =PS ( II) 

15  CONTINUE 

C 

RETURN 

25  WPITEC  6s  800 ) 

STOP 

C 

C FORMAT  STATEMENTS 

C 

800  FORMAT (///10X»20( 1 H* )///10X*  2 2 HE  XCE  SSI VE  ITERATIONS  / 

+ 1 OXo  9HIN  SHAFTP  /// 1 OX . 20 ( 1 H* )////// ) 

8 02  FOFMAT< //I IXo lHIo2X„2HIT# 1 2 X » 2HF I « I 3X.2HFN, 13X, 2HFP  « 1 2X  * 3 HD PI . 

+ 1 3X  o 2HDP* 1 2X  « 3HOPP  « 13X, 2HPI ,3X.5HPHASE  /) 

804  FOP MAT <8X. 2I4#3E15.4,4F15*6tI5) 

806  FORMAT ( ///10X*  6 H SHAFTP ) 

END 

tDHDGsP  * SUBROUTINE  PZAD  .L«l 


D25 


SUBROUTINE  PZAD 


l 


* 


aNBS*PLIBS.SHOW 

SUBROUTINE 


A. PZAD 
PZAD 


C 

C 

c 


2 


1 


8 

10 


1 5 
20 

800 

801 

802 

803 


THIS  POUTINE  CORRECTS  PZ  TERMS  FOR  PRESSURE 


PARAMETER  MM=220 . MS=2 , MC=4 . MP0=2 . MTP=2 * MFL=1 05 . MB= 1 05 
COMMON  NT,  P( MM ) «C( MM»MC) »NC( MM) • JC( MM*  MC) * ITS(MS1» 

1 FCtMMg  MC) « PZ (MMC MC)*PO(MM, MPO) tCOC  MM, MPO)  *F( MM ) # PFO( MFL# MPO) « 

2 FF ( MM ) e FO ( MM , MPO ) ,CS( MS)  ,PS<  MFL ) ,NS1 ( MS)  , NS2 ( MS ) • 

3 FSS(MS),N*NS  »NPO , I CONV , E , I BUG, X ( MM , MC ) » XO  ( MM , MPO) , TIT SHI  MS ,3 ) , 

4 NH,H<  MFL)  ,IFLOOR( MM)  , T ( MTP , MFL ) , NFS1 C MS) , NFS2<MS)t  IT (MB) ®NTP 

5 pNCO(MM)# JOCC MM, MPO) # TOUT 
COMMON  /PZZ/  PGZ 

DOUBLE  PRECISION  PgPOpPS 
I F ( IBUG  .GT.  — 2 ) GO  TO  1 
WRITE( 6, 800) 

DO  2 1 = 1 p N 
NN=NC .(  I ) 

I F ( NN  .EQ.  0 ) GO  TO  2 
W R I TE ( 6 , 80 1 )(  I , J , P Z ( I • J ) • J=1 ,NN ) 

CONTINUE 


NP 1 =N+ 1 

WRITE <6© 80 2) ( IL*PZ( IL* 1 ), IL=NP1 , NT ) 

DO  10  1=1 c N 
NN=NC< I ) 

I F ( NN  • EQ*  0 ) GO  TO  10 
I A= I T ( I ) 

I F 1= IFLOOR ( I ) 

DO  8 J J=1 » NN 
J=JC( Ip JJ) 

I FJ  = IFLOOR ( J ) 

I F( I F I • EQ • I F J ) GO  TO  8 
I B= IT ( J ) 

TEMPA=0.5*(T(  IA»IFI)+T(IB»IFJ)  ) 

PAVE=0.5*(P(I )+P( J) )+PGZ 

PZ(  I. JJ)=(0.03416*PAVE/TEMPA)*(H< IFJ)-H( IF  I ) ) 

CONTINUE 

CONTINUE 

DO  20  I S=1 , NS 

N 1=NS 1 ( IS) 

N2=NS2( I S ) — 1 
I TT = ITS ( IS) 

DO  15  I =N 1 , N2 
I FI= IFLOOR ( I ) 

IFJ=IF 1+1 

TEMPA=0.5*(T( ITT, I FI )+T< ITT.IFJ ) ) 

J=  I + 1 

PA  = 0. 5*  <P< I >+P<  J) ) +PGZ 

PZ(  I,  1 ) = <0.03416*PA/TEMPA)*(H(IFJ)-H<  IF  I)  ) 

CONTINUE 

RETURN 

p OP MAT (/10X, 10HINITIAL  PZ  /) 

FORMAT ( lOXp 3HPZ (, 12, 1H, 12, 4H)  = ,F12.4) 

p OR  MAT( lOXp  3HPZ ( , I 2 , 6H . 1 ) • ,F12.4) 

FORMAT ( / 10X , 1 1 HADJUSTED  PZ  /) 

END 
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5>NBS*PLIB$.SH0W  A.  OUT 
SUBROUTINE  OUT 
C 

c 

C THIS  ROUTINE  OUTPUTS  FLOWSAND  DIFFERENTIAL  PRESSURES 

C FOR  ALL  SHAFTS  AND  BUILDING  COMPARTMENTS 

C 

c 

PARAMETER  MM=220 * MS=2 . MC =4 ,MP0=2 « MTP=2 , MFL  = 1 05 . MB= 1 05 
COMMON  /CORR/C1  (MM , MC I , C2 < MM , MC ) • CO 1 ( MM, M PO) . C02 ( MM.  MPO ) 

COMMON  /IO/T ITLE( 12 ). I OUT. IUN I T. NC OMP ( MFL ) . SNC  OMP ( MFL ) 

COMMON  NT*  P(MM) . C(MM« MC) sNC(MM).  JC( MM.  MC) . ITS(  MS). 

1 FC  < MM. MC ) .PZ(MM. MC) . PO( MM. MPO) .CO(MM.MPO)  ,F( MM ) • PFO( MFL. MPO) . 

2 FF (MM ) , FO ( MM  , MPO ) *CS( MS)  .PS(MFL)»NS1 ( MS)  » NS2 ( MS ) . 

3 FSS(MS).N.NS  »NPO .ICONV.E.I BUG* X(MM»MC) »XO( MM  « M PO )»TITSH(MS»3) . 

4 NH,H( MFL ) , IFLOOR ( MM) ,T(  M TP, MFL ) .NFS1 ( MS)  . NFS2<  MS ) . IT ( M 8) . NTP 

5 cNCO(MM). JOC(MM, MPO). TOUT 
DOUBLE  PRECISION  P.PO.PS 
INTEGER  COM 

C 

C I UN  IT  = 1 FOR  SI  UNITS 

C I UN I T = 2 FOR  ENG  UNITS 

C WHEN  IUN  I T = 2 GO  TO  100 

I F ( I UNIT  .EQ.  2 ) GO  TO  100 

C 

C BUILDING  COMPARTMENT  OUTPUT 

C 

1  1=0 
I L=  0 

WRITE (6.800) ( T I TLE ( I)  .1=1  .12) 

DO  30  IFF=1,NH 
NNN=NCOMP( IFF ) 

I F ( NNN  .EQ.  0 ) GO  TO  30 
DO  29  I C=1 . NNN 
1 = 1 + 1 
NN=NC ( I ) 

NNC=NCO( I ) 

I L=  IL+NN+NNO+2 

I F ( IL  *LT.  51  ) GO  TO  2 

WRITE( 6. 800 ) (TITLE!  I ) .1  = 1 .12) 

I L=  NN  + NNO+  2 

2 IF(NN  .GT.  0)GO  TO  3 

WR ITE(6.801 ) IFF ,IC,P( I ) , IT( E ) ,FF(I ) 

GO  TO  21 

3 DO  20  J=1 . NN 
J J = JC ( I. J ) 

DP  = P ( J J ) — P ( I ) + PZ< I , J) 

CC=C2 ( I. J ) 

I F ( DP  .LT.  0. ) CC=C 1 ( I . J ) 

I F ( J J .LE.  N ) GO  TO  10 
DO  5 I S= 1 . NS 

I F ( J J .GE.  NSl(IS)  .AND.  JJ  .LE.  NS2(IS))G0  TO  6 

5 CONTINUE 

6 IF(J  .GT.  l)GO  TO  7 

WR ITE(6,802) IFF,IC.P( I ) , I T( I ) ,FF( I ) . (TlTSH(  IS.K  ) . K= 1 .3) 

+ ,DP.CC«X( 1,1  ) * FC (1,1) 

GO  TO  20 
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* 

7 WRITE (69 80  3) ( TI TSH ( I S* K ) • K= 1 1 3 ) • DP.  CC • X ( I . J ) . FC ( l • J ) 

GO  TO  20 

10  IFJ=IFLOOP( JJ) 

COM=J J-SNCOMP ( IFJ ) 

I F ( J • GT • 1 ) GO  TO  12 

WRITE(6,804)IFF»IC»P(I )• IT( I) »FF ( I ) . I F J ,C OM«DP , CC , X ( I . 1 ) ,FC( I .1 ) 

GO  TO  20 

12  WRITE (6. 805) I F J .COM , DP , CC  tX ( I , J) ,FC( I. J) 

20  CONTINUE 

21  I F ( NNO  • EQ • 0 ) GO  TC  29 

DO  23  J=1 « NNO 
JJ=JCC< It J) 

DP=PO(  Ic J )-P(  I ) 

CC  = C02< It  J) 

I F ( OP  .LT.  O. )CC=C01  ( I • J ) 

23  WRITE(6»806)JJ»  DP  ,CC,XO(I  «J)«F0(I,J) 

29  WRITE(6# 807)F( I ) 

30  CONTINUE 
C 

C SHAFT  OUTPUT 

C 

IL=  IL+6 
DO  60  IS=1,NS 
N1=NS1 ( IS) 

N2=NS2( IS) 

IL=  IL+N2+ 1 8— N 1 
I F ( IL  .LT • 56 ) GO  TO  32 
WP I TE (6*814) 

IL=N2+18-N1 

32  WRITE (6, 308) (TITSH( IS.K). K=1 .3)«ITS(IS).CS(IS) 

DO  50  I=N1,N2 
NN=NC( I ) 

IF ( NN  .GT,  O ) GO  TO  35 

WRI TE( 6,809) I FLOOR ( I ) * P ( I ) , FF ( I ) 

GO  TO  41 

35  DO  40  J = 1 9 NN 
JJ=JC( Ie J) 

DP=P( J J )— P ( I ) 

CC=C2( 1 9 J ) 

I F ( DP  .LT.  0. )CC=C1(I,J) 

I F J = I FLOOR ( J J ) 

COM= J J— SNCOMP ( IFJ ) 

I F ( J .GT.  1 ) GO  TO  36 

WRITE(6»810) I FLOOR ( I ) , P ( I ) , FF ( I ) , I F J , COM, D P ,CC , X ( I , 1 ) , FC ( I , 1 ) 

GO  TO  40 

36  WRI  TE (6o 81 1 ) IF J,COM,DP,CC  , X ( I , J ) , F C ( I , J ) 

40  CONTINUE 

41  NNO  =NCO ( I ) 

I F ( NNO  .EQ.  0 ) GO  TO  50 
DO  46  J = 1 , NNO 
JJ= JCC( I. J ) 

DP=PO ( I c J >-P(  I ) 

CC=C02 ( I . J ) 

I F ( CP  .LT.  0. )CC=C01 ( I , J ) 

4 6 WRITF(6,812) J J , DP , CC , XO ( I ,J)  ,FO(I, J) 

50  CONTINUE 

WRITF(6c813)FSS( IS) 

60  CONTINUE 
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GO  TO  165 

c 

C BUILDING  DATA  OUTPUT  FOR  IUNIT  = 2 

C 

100  1=0 
I L = 0 

WRITE! 6,800) (TITLE! I) ,1=1 ,12) 

DO  130  IFF  = 1 « NH 
NNN=NCOMP< IFF ) 

IF! NNN  • EQ • 0 ) GO  TO  130 
DO  129  I C= 1 » NNN 
1=1  + 1 

FF I =F ( I )/0.47l9 
PI  I I = P( I ) /248 • 8 
FFF  =FF ( I J/0.4719 
NN=NC ( I ) 

NNC=NCO ( I ) 

I L=  IL+NN+NNO+  2 
I F ( IL  ,LT • 51  ) GO  TO  102 
WRITE (6,800) (TITLE! I)  ,1=1  ,12) 

I L = NN+NNO  + 2 

102  IF! NN  • GT • 0 ) GO  TO  103 

WP ITE(6, 60 1 ) I FF , I C , P 1 1 I , I T(  I) ,FFF 
GO  TO  121 

103  DO  120  J=1 ,NN 

F CC C=FC ( I e J ) /O • 47 1 9 
JJ=JC(  I, J) 

0P= ( P ( J J ) — P( I )+PZ( I, J))/248.8 
CC=C2 ( I* J ) 

IF! DP  ,LT.  0, ) CC=C 1 ( I , J ) 

CC=CC*2, 119/(0«004019**X( I , J ) ) 

I F ( J J ,LE.  N ) GO  TO  110 
DO  105  IS=1 ,NS 

IF! J J .GE,  NSl(IS)  .AND.  JJ  • LE • NS2(IS))G0  TO  106 

105  CONTINUE 

106  IF! J • GT • 1 ) GO  TO  107 

WRITE(6,602) I FF , IC,PIII ,IT( I ) , FFF  ,(TITSH( IS,K) ,K=1 ,3) 

+ »DP,CC*X(I ,1 ) ,FCCC 
GO  TO  120 

107  WRI TE! 6*603) (TIT SH ( I S®  K ) , K= 1 , 3 ) ® DP , CC • X ( I , J ) , FC CC 

GO  TO  120 

110  IFJ=TFLOOR( JJ) 

COM=J J—SNCOMP ( IFJ ) 

IF! J • GT • 1 ) GO  TO  112 

W R ITE! 6, 60  4)IFF,IC»PIII,IT(  I), FFF  ,IFJ,C0M,DP,CC,X(I,1),FCCC 
GO  TO  120 

112  WPITE(6»  605 ) IFJ, COM, DP, CC  »X(I,J),FCCC 

120  CONTINUE 

121  IF! NNO  • EQ • 0 ) GO  TO  129 
DO  123  J=1 ,NNO 
FCO=FO(  I, J ) /O ,471 9 

J J= JCC! I • J ) 

0°= (PO(  I« J)-P(  I ) )/248,8 
CC  = C02( I®  J ) 

I F ( CP  ,LT,  0 • ) CC=C01 ( I , J ) 

CC=CC*2. 119/! 0,  0 04019*^X0(1  ,J)  ) 

123  WRITE!6,606)JJ,DP,CC.X0(I,J) ,F00 
129  WRITE(6»807)FFI 
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130  CONTINUE 
C 

C SHAFT  OUTPUT  FOR  IUNIT  = 2 

C 

IL  = IL  + 6 

00  160  I S=  1 , NS 
CSS=CS( IS) /O • 02992 
FF I=FSS( IS) /O  *4719 
N 1 = N S 1 C IS) 

N2=NS2< IS) 

IL= IL+N2+18-N1 
I F ( I L .LT.  56 ) GO  TO  132 
WP ITE<6,81 4) 

I L = N2+ 1 8— N 1 

132  WRITE  (6s>808)<  TITSHI  ISBK)#K=1*3)«ITS(IS),CSS 
DO  150  I =N 1 « N2 
FFF  =FF ( I ) /O • 4719 
PII  I = P(  I)/248.8 
NN=  NC ( I ) 

I F ( NN  » GT • 0 ) GC  TO  135 

WRITE ( 6*609) I FLOOR! I ) *°I 1 I * FFF 

GO  TO  141 

135  DO  140  J=1 iNN 
FCCC=FC( Is  J)/0.4719 
J J = JC ( I e J ) 

DP=(P( JJ )-P( I ) )/248.3 
CC=C2 C I«  J) 

I F ( DP  .LT,  0 • ) CC=C 1 ( I * J ) 

CC=CC*2. 1 19/< 0. 00401 9**X( I, J) ) 

I FJ= I FLOOR ( J J ) 

COM= J J— SNCGMP ( I F J ) 

I F ( J .GT.  1 ) GO  TO  136 

WDITE(6»61 0) IFL  00 R ( I ) » P 1 1 I » FFF  , I F J , COM, DP , CC  , X ( I , 1 ) , FCCC 
GO  TO  140 

136  WRI TE(6c61 1 ) IFJcCOMo  DP  c CC  « X(I*J),FCCC 

140  CONTINUE 

1 41  NNC  = NCO(  I ) 

I F ( NNO  .EQ.  0 ) GO  TO  150 
DO  146  J= 1 , NNO 
F 00  =F0 ( I, J )/0.4719 
J J=  JOC ( Is J ) 

DP=(PO(I*J )— P ( I ) ) /248. 8 
CC  = C02 ( Io  J ) 

I F ( DP  .LT.  0. ) CC=C01 ( I , J) 

CC  = CC*2.  1 19/(  0. 004019**XO( I • J)  ) 

146  WPITE(6»612)JJ»DP»CC»X0(I »J) ,F00 
150  CONTINUE 

WRITE(6«  813) FF I 
160  CONTINUE  < 

C 

C SUMMARY  OUTPUT 

C USER  INSERTS  WRITE  STATEMENTS  TO  FILE  I OUT 

C 

165  CONTINUE 
RETURN 

C 

c 

C FORMAT  STATEMENTS 
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C 

60  1 
602 

603 

604 


60  5 
606 
609 
6 10 

611 

612 

800 


FORMAT !/4X, I3,I10.F13.3,I8.F12,0)  STURHI .A  W0H8 .$ P I JR*8 GUG 

FOR  MAT  ( /4X  * I3,110,F13.3,I8,F12.0,3X  ,3  A6  , i F ?8iI)VP  f 8^  £ »F 1 1 • 1 ) 

FORM AT { 53X, 3 A6 , FI  6. 3, FI  5.0, FI  0.2, FI  1. 1 ) 

FC5GH»?3/3*«f  3#f  fe&M^ARTMENT*  13 • 

i fi  i.3afT6tio9?i0J€9Pi¥Iiy  obriudbr  si  hur  rbhtoha  31  ro 

FOR  MAT  ( 53X»  5HFL00R  « 13*  12H  COMPARTMENTS3?^**!'  1 ?§  .Pfd.tyjP  10  ,2  »F1  ? .1  ) 
F OR  MAT  ( 53  X » 1 7 H OUTS  I DE  D IRECT  ION  • 13  , FI  4 . 3#  F ?3  .80^1  6 • 2#  1 • 1 ) 

FORMAT < 4X® I 3»  F 1 0.  3 *F  1 1 • 0 ) 

FORMAT!  4X,  I 3 . F 1 0 .3  , F 1 1 • 0 , 3X  . 5HFL00R  , 1 3 • 1 aH^Otf^RfMEN?'?  13  ,F11 .3, 

1 F15.0.F10.2.FU.1  ) YJMQ  2THI0R  S H33WT38  3 ROB  4 = D» 

F0RMAT(31Xt  5HFLOOR  *13*1 2H  COMPARTMENT, I3,F11.3.F15.0,F10.2,Fll2l) 

fop m At  r^5MPl?Xg§*2  8x« 

* IGR*  •43HMi*RH*eWt¥A49SH*«Mr*QawtemiWe^4R¥9MNL§0Rl2H»¥tfi>«gaM*4RTMENT 

2 ,2x.8HP*eBfcdfte«2Ritf?APMcMettua«4fti*&e9f s,  To  , 

. <c*23*eRT«BPAt9!Swms;9*#itHcHB^iWei^T?aXV§fl£xP0NEfl¥#‘SxUMi>&0W  S') 

8 0pTM^dR»i*tt^4Ma$3?¥l<jSFi3Sf^lfeVP¥29^>T,  CMM  >fl00J3I.  ( J3M  >H,HH  4 
80  2 FORMAT  !/4X»  13, 110, F 13.1,  1 8 » F 1 2,  <6  V§*  9 i ?<?£  i 2 SF 1 1 .1  ) 

F0RMAT(53X,3A6,F16.1 ,F15.1 , FI 0 .2 * I Of » q H0ISI03RR  3J0UOO 
FOi*Ma¥IS»¥«?6*3#  ^9M^»6tment#I3, 

1 FI  1. Is F15.1* F10.2«F11 . 1 ) HURI\HUR\  HDMMOD 

FORMAT! 53XC5HFLOOR, I 3, 12H  COMPARTMENT,  I 3,  FI  1 . 1 tPVS  i Pf  b .2  ,F1  1 .1  ) 
FORMAT (53X* 17H0UTSI0E  DIPECTI0N,I3,F14.1,F15.1, FI 0 . ¥S  i 1 ) 

FORMAT! 1 15X,F8,1 .4H  NET)  MHW( OOY «5 ) 0 A3 R 

FOFMAT(////20X,  3A6//20X  ,2  OHTEMPERATURE  !£?  • T v>  i 3^2  d H i 

1 238-SHAFT  FLOW  COEFF I Cl  0N?I  • 1 £ f 9 *?8W  AB.il&SfgOi*  £3*^  SNF I XED  , 

2 2 8 X ® 1 2 HD  I FFERENT  I AL  » 5 X*  4 HFLOW,  8Xe  4HFLQW/4  WSHpEdSR&g*,  9GPRESSURE, 

3 5X«  4 HFLOW  «3X«  1 6HCCNNECT  I ON  TO , 1 2X  , 88Pfc£s3tf fcEQSx S i ?At0 EFF I C I EN  T 

4, 2X ,8HEXP0NENT,5X,4HPQ§w0?)0D! 4 .TO,  OX  .RO . I .TJ.  0X)3I 

FORMAT  !4X,  13, F10, 1 , FI  1 .0)  3=88 

FORMAT!  4XS  I 3,  FI  0,  1 , F 1 1 • 0 • 3X  , 5HFL00R  , f>  3 , 91  H°§<*>tfp  a6¥  MElQ¥  ) ? 3 , FI  1.1, 

1 FIE.  1 .F10.2.F1 1 .1  ) °S  OT  00 ! S .03.  OX)BI 

FORMAT  !31  Xs  5HFLC0R*  13  s 12H  C OMPA  RTME  Nt  f I SjrF??^  fc  F?§  • 1 id  J 2,  F 1 1.1) 

F OP  MAT  !31X®  1 7H0LT  SIDE  DIRECTION  s I 3 , P?  4S?f  , P?  4$  1 tPl  © • 2$  I*  I i • 1 ) 

FOR  MAT! 93X„F8. 1 ,4H  NET)  ROTS 

FORMAT! 1 HI  ) 

END  21  T3AHS  HI  33  HI  30HAHD 


803 

804 

805 
8 0 6 

807 

808 


909 

810 

81  1 
812 

813 

814 
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SUBROUTINE  IN°UT  £2 I ) Igfclf! H 

(SI )£3H=SH 
P IY4. 0*88=08 ( S .03.  TIHUI)3I 
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SUBROUTINE  INPUT2 


* 


£>N3S*PL  IB*.  SHOW  A.INPUT2  « 0 • S I B , 8 I * F . 0 I -I  , Z>  i I » " 

< I . I : SUBROUTINE  Xygyi^.dAE.  XE»  0 .£  13 , B I • C *E  13  *0 1 I*  £ I *XAN  ) TAM 

c (X .1 13, S. 013,0.6 I3,0.dI3.dAE,X06> TAMPO^  EOd 

C THIS  ROUTINE  DETERMINES  IF  THE  PROGRAM • IS  TO  EE  STOOGED 

C OR  IF  ANOTHER  RUN  IS  REQUIRED  WITH  MODIFIED  DATA 

C ilWW.  £9*1*  if  ifN'MJf e?th3MTPAPM00  HS  f *EI  . POO  JBH6  *XE6  ) TAM  < 

c r * I «E.  A 13  ,01  .MOITDBPIC  BO  I 3TU0HT  I * XE8) TAMPQ3 

C KC  = 2 FOR  C ( 0 • 1 1 3 » 0 • 0 I 3 o F I » • '. 

,€.1X3,  0I$fH5M?p£9Bo§oHSI .0  I*  P00J3HB . XO, 0.113,0.0X3, 01  »XA  ) T A M 

C KC  = 4 FOR  C BETWEEN  2 POINTS  ONLY  ( I . I IB  ,S*  0! 

<lClI3,S.0I3,0.6I3,0.1I3.0I .TMBMTBAPMOD  HS 1 ,0  I , BOO J3H6 * X I E > T AM  303  X Id 

PARAMETER  MM=220 , MS=3f».*IC=*  ;.MPE)=a»  MTReSsMfrfeR*  ©S«M8b  106 
COMMON  NT*  P(MM|x§fMMfM^,NE<;MMi#UG(MM*MC>f  fTS(MS)* 

1 FC(MMoMC>*  PZ<MM#  MC*»P9(;MM»MPOJf§9ftMM«MRQiAf  < MM  >*  PFO<  MFL*  MPO)  • 

2 FF  (MM)oFO<MM,MPOIfCS(iMS)  ,R§  ( MRMa  NSli  t MS*  t MSB  <MS»  * 

3 FSS<MS)»NoNS  *NPO* ICON V*  E * I BUG  > X ( MM  «MC ) *X  0( MM,MP0> .TITSHCMS.3) . 

A NH  * H(  MFL  ) . IFLOOR  ( MM ) * T < f f MSk  I W?8?  t MS  ) ♦ IT  ( M B > . NTR  0 8 

5 * NCG<  MM),  JOCCMM*  MPP.£t*eM¥.S  13,81  , I „EIB>OX  I : ’ I v 

DOUBLE  PRECISION  P « ?0 f»?$ , S . 0 13 , I. 3X3, I • d 13 • dAE# XEB> TAMPO B 
COMMON  /IO/THiyE(Jta^»XeW¥»l9NI¥#NC0MR€MFW) *SNCOMP(MFL) 

COMMON  /RUN/I  RUN  (IeII3»S»OI3*I.6I3oX.I  IB 

X.  X XB, s. 89M9?169Ntl . I ib,oi  , tmbmtpapmod  HS  f ,0  I ,RQOJ3HS «X05 itampo 
I lf¥^T§*OIB ,I.5I3,i.AI3,0I,M0ITD3BIG  HCI6TU0HYI  «XE6)TAM  303 
5 READC  5*700  )NNN  ( T3H  H A » X • 8B  » X 6 X X ) T A 

lRtHN«iE«9T*B9!850R  3PUTAP3qM3TM0S,  XOSWdAO  , XOSWW  ) TAM  BOB 
WRI TE<6, 600) (TITLE! »lBl *12)  • I RUN  D IBB  BOD  WCJB  / -i  “ 

IU8  2BPBR8  **§Qb$&3%WB9Wa\W0J3HA  »X8  « WO  JBH  A»X3#JAI  TH3B333  I CHS  I BX  8< 

T MB I D IBBBPF AC (5c7OO)KC»IS*08,XSX,OT  HOI TDBMHOOHd  X *XE«W0J3HA » X?  0 

I F ( KC  »LT • 1 • OR*  KC  • GT • 4 ) GO  ( ^OwgftOlHA  , XE  , TM3H0PX3H8  , XS  » * 


BB=e  (0* X XB» X .0X3 „0I  ,XA ) Tt 

, X . X X^tfcft  ^ft3MfflARfcd§0Hl9,l;Q.  ROOJ3HE,XO.O  • 113,1.013,01  ,XA  ) TAMPO  B 
I F{  KC  .EQ.  2 ) GO  TO  20  ( X . i X B , S • 0 XB , X . 3 X B X 

X • X X3,sl6f*tr«§9R.n§93't@l39MBMTPAqMOO  HS  X *01 »ROO  JBHE  8X IE) TAMPO - 
• If  I 9 HOI  T 03  3 1 Q BO  I 2 T JOHY  I a X IE)  TAM  3 

STOP  (T3H  HA, X • 83  * XO 9 ) TA M 3Q3 

( XHI ) T AMPOB  AIR 

CHANGE  IN  FF  IN  SHAFT  IS 
0 N1=NS1(IS)  3HITU0P3US 


N2=NS2( IS) 

I F ( I UN I T .EQ.  2 ) BB=BB*  0.4719 
DO  15  I=N1,N2 
FF ( I )=BB 
5 CONTINUE 

WRITE! 6c  800) I S, B 
GO  TO  100 

CHANGE  IN  C FOR  SHAFT  IS 
C 

20  N 1 =NS 1 ( IS ) 

N2=NS2< IS) 

DO  28  I=N1,N2 

NN  = NC ( I ) 

I F ( NN  .EQ.  0 ) GO  TO  28 
DO  24  J=1,NN 

IFCIUNIT  .EQ.  2 ) BB  = B*0.47  19/( 2A8.8**X( I • 1 ) ) 

C ( I c J )=BB 
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INPUT2 


* SUBROUTINE 

JJ=JC«I*J) 

NM=NC( JJ) 

30  22  K = 1 e NM 

IF<  JC< JJ.K ) .EQ.  I ) GO  TO  23 

22  CONTINUE 

23  CIJJ*K)=8B 

24  CONTINUE 

26  CONTINUE 

WRITE<6* 80i)IS,B 
GO  TO  100 

c 

C CHANGE  CO  IN  BUILDING 

C 

30  IF  I IS  • GT  * 0 ) GO  TO  40 

30  38  1=1. N 
NNO=NCO ( I ) 

I F < NNG  • EQ.  0 ) GO  TO  38 

00  34  J= 1 . NNO 
IF«IUNIT  ,EQ.  2)BB=B*0 .4719/(248.8* *XO<I»J)) 
COd.J  ) = BB 

34  CONTINUE 

38  CONTINUE 

WRITE(6e802)B 
GO  TO  100 
C 

C CHANGE  CO  IN  SHAFT  IS 

C 

40  N 1 = N S 1 ( IS) 

N2=NS2( IS) 

DO  48  I=N 1 • N2 
NNO  =NCO ( I ) 

1 F ( NNO  .EQ.  O ) GO  TO  48 
00  44  J = 1 » NNO 

IF*  IUNIT  .EQ.  2 )BB=B*0.4719/(248.8**X0< I, J)) 

C0{  Ic J)=BB 
44  CONTINUE 

48  CONTINUE 

WP  I TE ( 6 ® 8 03 ) IS.B 

GO  TO  100 
C 

C NEW  C BETWEEN  POINTS  IS  AND  JS 

C 

50  PFAD(5»700)JS 
NN=NC* IS) 

DO  51  J = 1 e NN 

IF*  JC ( IS. J ) .EQ.  J S ) GO  TO  52 

51  CONTINUE 
GO  TO  200 

52  IF(  IUNIT  .EQ.  2 ) BB=B*0 • 47 19/< 248 • 8* *X(  I S . J ) ) 
C( I Sc J )=BB 

NN=NC( JS) 

DO  53  J= 1 c NN 

I F ( JC ( JSc  J ) .EQ.  IS) GO  TO  54 

53  CONTINUE 
GO  TO  200 

54  C(JS«J)=BB 

WRITE  < 6c  805)  IS, JS, 8 
100  CONTINUE 
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200 

600 

700 

800 

801 

802 

803 

804 

805 


RETURN 

WR I TE (6#  804 ) 

STOP 

FORMAT! ///I OX.  12A6//10X. 3 HR UN. I 4/// ) 
F OF  MAT ( ) 


FORMAT! 10X.25HNEW 
FORMAT! 1 OX, 16HNEW 
FORMAT! 10X.23HNEW 
F OF  MAT! 10X.  17HNEW 


FIXED  FLOW  FOR  SHAFT  ,I2,4H 
C FOR  SHAFT  *12, 4H  IS.F8.1//) 
CO  FOR  BUILDING  IS  »F8*1//) 

CO  FOR  SHAFT  .I2.4H  IS.F8.1// 

F0RMAT!//10X. 2 1 HERR OR  IN  INPUT 2 DATA  /) 

FORMAT! 10X. 14HNEW  C BTTWEEN  ,I3.4H  AND.I3.3H  IS 
END 


3)  HD  Go  ° * 


SUBROUTINE  UNITS 


IS.FS. 1/) 
) 

,F8.1//> 

• L , 1 
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SUBROUTINE  UNITS 


* 


S)NBS*PL  IBS.  SHOW  A. UNITS 
SUBROUTINE  UNITS 
C 

C THIS  ROUTINE  CONVERTS  VARIABLES  H,FF,C,CO,CS  TO  SI  UNITS 

C 

c 

PARAMETER  MM  = 220 . MS=2 « MC=4 , MP0=2 • MTP  = 2 • M FL=1 0 5 . MB= 1 0 5 
COMMON  NT,  P( MM) ,C( MM,MC) ,NC( MM) . JC( MM. MC) , ITS( MS ), 

1 FC<  MM e MC) ePZ(MM» MC ) , PO  ( MM» MPO ) * CO ( MM. MPO) *F( MM). PFO( MFL.  MPO) . 

2 FF (MM) ,FO< MM, MPO) ,CS( MS) , PS ( MFL ) ,NS1 ( MS) ,NS2(MS) , 

3 FSS(MS).  N.  NS  ,NPO,  I CON  V , E , I BUG  , X ( MM  , MC  ) , X O ( MM  , M PO  ) , T I T SH  ( MS  , 3 ) . 
A NH  *H ( MFL ) , I FLOOR (MM)  . T ( MTP , MFL ) . NF SI (MS)  ,NFS2(MS),  IT(MB).NTP 

5 . NCO(MM ). JOC( MM, MPG) • TOUT 
DOUELE  PRECISION  P,PO,PS 
DIMENSION  B ( 4 ) 

DATA  B/0* 30 48 ,248. 8,0.4719, 0. 02992/ 

DO  10  1 = 1.  NH 

10  H( I )=H( I )*B( 1 ) 

DO  20  1=1, NT 

FF ( I )=FF( I )*B(3) 

DO  16  J= 1 , MC 

IF(X(IeJ)  .EQ.  0 ) GO  TO  16 
BB=E (3  )/(B(2 )**X(  I, J)  ) 

C(  I * J ) =C ( I, J)*BB 
16  CONTINUE 

DO  18  J= 1 « MPO 
IF(XO(I®J)  .EQ.  0 ) GO  TO  18 
BB=B(3)/(3<2)**X0( I , J) ) 

CC(  I g J ) =C0 ( I , J )*BB 
18  CONTINUE 

20  CCNTINUE 

DO  22  I S= 1 . NS 
22  CS(  IS)  = CS(IS)*B(4) 

RETURN 
END 

* 


£HDGe  =» 


SUBROUTINE  SIMEQ 
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SUBROUTINE  SIMEQ 


* 


©NB  S*PL I B $ • SHOW  A*  SI MEQ 

SUBROUTINE  SIMEQ 
C 

C CHOLESKY*  S METHOD  OF  SOLUTION  OF 

C SIMULTANEOUS  LINEAR  ALGEBRIC  EQUATIONS 

C 

PARAMETER  MM=220, MS=2*  MC=4* MPO=2 • MTP  = 2 ,MFL=1 05 • M8=l 05 
PARAMETER  MBP=MB+ 1 
DOUBLE  PRECISION  A*X 
COMMON  /MAT/  A< MB • MBP ) , X< MB ) « N 
NP1=N+1 
ZEP.G=  1 • 0E—35 
K = Q 
C 

c 

C SEE  IF  A ( 1 • 1 ) IS  ZERO 

C IF  SO  ADD  ANOTHER  ROW  TO  ROW  1 

IF( ABS( A( 1 ♦ 1 ) ) .GT.  ZERO) GO  TO  40 
DO  31  1=1, N 

1 F( A( I« 1 ) • NE • 0* ) GO  TO  32 

31  CONTINUE 

12  WRITE( 6, 804)K 
STOP 

32  DO  33  J= 1 « NP 1 

33  A ( 1 o J )=A< 1, J ) +A ( I, J ) 

C 

C CALCULATE  UPPER  AND  LOWER 

C TRIANGULAR  MATRICES  OVER  ORIG 

C MATRIX  A 

40  A A=  A ( 1 c 1 ) 

DO  2 J=2, NP1 

2 M1«J)=A(1,J)/AA 
DO  10  1 = 2*  N 

K = 0 
C 

C STORE  A ( I o 1 ) •••  A ( I , I ) IN  X ARRAY 

C IN  CASE  NEW  A(I,I)  IS  ZERO 

C ROW  I CAN  BE  RECALCULATED 

4 DO  F J = 1 * I 

5 X ( J ) = A ( I * J ) 

K = K+  1 

DO  10  J=2*NP1 
I F( J • GT • I ) GO  TO  8 
J Ml  = J- 1 
A 4=0. 

DO  3 IR=  1*  JM1 

3 AA=AA  + A(  I*  IR) * A ( I R , J ) 

A ( I * J ) = A ( I, J )-AA 

C 

C CHECK  IF  A ( I , I ) IS  ZERO 

C IF  SO  MULTIPLY  OLD  ROW  I BY  2* 

C 

I F < I .NE,  J ) GO  TO  1 0 

IF( ABS(A(I * I)  ) .GT.  ZERO ) GO  TO  10 

DO  6 JJ=1,I 

6 A ( I * J J ) =X ( J J ) 

DO  7 J J= 1 * NP 1 
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C 

c 

c 

c 

r 


c 


7 M I o J )=2.*A( I,  J ) 

I F ( K • GT • 3 ) GO  TO  12 

GO  TO  4 

8 I M 1 = I—  1 
A A — 0 * 

DO  9 I R=  1 * 1 M 1 

9 AA  = AA  + AU«  XR)*AC  IR»  J) 

A (I#  J)=( A(I ,J)-AA)/A(1 , I) 

10  CONTINUE 

END  OF  CALCULATION  OF  TRIANGULAR  MATRICES 


BACKWARD  SUBSTITUTION 

X(N  ) = A(N*NP1 ) 

DO  20  I 1=2 • N 
A A=0« 

I =NP1-I I 
IP1=I+1 
DO  15  J=IPl,N 
15  AA  = AA  + A<  I®  J)*X( -J) 

20  X ( I ) = A(  I • NP1 ) — A A 

804  ^ OF.  MAT(  //////l  OX  » 1 6HPR0GR  AM  FAILURE  ,13/////) 
END 


fO  HD  Go  P * 


FUNCTION  FLOW  «L,1 
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FUNCTION  FLOW 


* 


5>NB  S*°L  I B $•  S HOW  A. FLOW 

FUNCTION  FLOWCPI  ,PJ,PZ.C,X) 

OOUELE  PRECISION  PI.PJ 

THIS  FUNCTION  CALCULATES  FLOWS  BETWEEN  TWO  POINTS 

C 

I F ( ABS(C)  »LT  • OtOODGO  TO  10 
C>o=p  j_pi+pz 
S I GN= 1*0 

I F ( OP  .LT.  .0)SIGN=-1. 

I F ( X .NE.  *5 ) GO  TO  11 
FLOW=SIGN*C*SQRT< SIGN* DP) 

RETURN 

10  FLOU=0.0 
RETURN 

11  FLOW=SIGN*C*< < SIGN*DP)**X> 

RETURN 

END 

SHDG.P  * FUNCTION  PFLOW 


• Lei 


D38 


ro 


FUNCTION  PFLOW 


3)NBS*PL  I 8$  • SHOW  A. PFLOW 

FUNCTION  PFLOW ( I * P I ) 

c 

c 

C THIS  FUNCTION  CALCULATES  NET  FLOWS  INTO  POINT  I 

C 

PARAMETER  MM=220 . MS=2. MC=4 ,MP0=2,MTP=2,MFL=105 ,MB=105 
COMMON  /CORR/C1 (MM.MC) ,C2 (MM * MC ) • COl ( MM* MPO) , C02 ( MM, MPO ) 

COMMON  NT o P ( MM ) ♦ C ( MM • MC) • NC( MM ) • JC( MM* MC ) ♦ ITS ( MS ) • 

1 FC(MMoMC) * PZ ( MM  * MC) »PO(MM»MPO) *CO( MM, MPO)  ,F ( MM ) , PFO( MF L * MPO)  , 

2 FFCMM) ,FO( MM.MPO) ,CS( MS) .PS(MFL) ,NS1 ( MS) ,NS2(MS) , 

3 FSS ( MS ) * N»  NS ,NPO, ICONV, E » I BUG ♦ X ( MM , MC ) ,X 0<  MM , M PO) , TI TSH<  MS , 3 ) , 
A NH , H ( MFL ) , IFLOOR(MM) , T ( MTP , MFL ) , NFS1 (MS) , NFS2(MS), IT(MB),NTP 

5 e NCO< MM) , JOCC MM, MPO) ,TOUT 
DOUBLE  PRECISION  P,PO,PS.PI 
NN=NC( I ) 

SUM=Q • 

IF(NN  «EQ.  0 ) GO  TO  3 
DO  1 JJ=1,NN 
J=JC< I, JJ ) 

CC=  C 1 ( I , J J ) 

IF(PI  ,LT,  P ( J ) ) CC=C2(  I » J J ) 

PZZ=PZ( I , J J) 

I F ( I • GT • N ) PZZ=0 • 

F C ( I , J J ) =FLO W (P  I « P ( J ) »PZZ»CC»X( I,JJ ) ) 

1 SUM=SUM+FC( I,  J J ) 

NNC=NCO( I ) 

I F ( NNO  • EQ  o 0 ) GO  TO  4 
DO  2 K= 1 , NNO 
CC=C01 ( I,K ) 

I F ( P I .LT.  PO (I  , K ) ) CC=C 02 ( I * K ) 

F0( IcK ) = FLOW(P I ,P0  U ,K ) ,0 ,CC, X0( I ,K) ) 

SUM  = SUM  + FO( I * K) 

PFLOW=SUM+FF( I ) 

I F ( I ,LE,  N)F( I )=SUM+FF( I ) 

RETURN 
END 

BHDG,P  * FUNCTION  PSFLOW  ,L,1 
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a)N9S4'PL  I B $• SHOW  A.PSFLOW 

FUNCTION  P SFLO W ( I » P I ) 

C 

c 

C THIS  FUNCTION  CALCULATES  NET  HORIZONTAL  FLOWS 

C INTO  A FLOOR  OF  A SHAFT 

C 

PARAMETER  MM=220 , MS=2 . MC=4 ,MP0=2, MTP=2 . MFL=1 05 .MB=105 
COMMON  /C0RR/C1  (MM , MC) * C2 (MM »MC ) • C01 C MM* MPO) • C02( MM, MPO ) 

COMMON  NT,  P(MM)*C( MM*  MC) * NC (MM), JC( MM*  MC),  ITS(MS), 

1 FC(MM,MC) *PZ ( MM,MC) ,PO( MM, MPO)  ,CO( MM, MPO)  , F ( MM  ) , PFO( MFL, MPO ) , 

2 FF(MM) e FO ( MM  * MPO ) ,CS( MS) *PS(MFL)*NS1 (MS) ,NS2(MS) • 

3 FSS(MS) *N*NS,NPO. I CON V , E , I BUG , X ( MM • MC ) • X 0 ( MM , M PO ) • TI T SH( MS ,3 ) , 
A NH ,H ( MFL ) , IFLOOR ( MM ) ,T ( MTP , MFL ) . NFS1 (MS) ,NFS2(MS), IT(MB)*NTP 

5 cNCO(MM)*  JOC ( MM*  MPO ) 

DOUELE  PRECISION  P,PO,PS,PI 
NN=NC( I ) 

SUM=0* 


I F ( NN  *EQ*  0 ) GO  TO  3 
DO  1 J J= 1 * NN 
J = JC(  Ie  JJ) 

CC  = C1 ( I*  JJ ) 

I F ( P I ,LT,  P ( J ) ) CC  =C2( I »JJ) 

1 SUM=SUM+FLOW(PI ,P ( J ) ,0 , CC ,X ( I • J J ) ) 

3 NNC  =NCC ( I ) 

I F ( NNO  • EQ*  0 ) GO  TO  A 
DO  2 K= 1 , NNO 
CC=C01 ( I,K ) 

I F ( P I ,LT*  P0( I ,K)  )CC=C02(I ,K) 

2 SUM=SUM+FLOW(PI ,P0 ( I , K ) , O , CC , XO ( I • K ) ) 

A P SFLO W=SUM+FF ( I ) 

RETURN 

END 


5>PRKPT  PRINTS 
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18.  Number  of  source  program  state- 
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